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INTRODUCTION 


CS" far as can be discovered, there is yet no account of times of 


ossification in the Australian aborigine. The present study, despite 


its obvious defects, thus gains some interest from that point of view. 

These observations are confined to the members of the Ngalia tribe, 
together with a very few pure or mixed offspring of neighboring tribes. 
They were all assembled as a single economic unit in and around the 
Government Settlement at Yuendumu, some 200 miles northwest of 
Alice Springs. The “aliens” in the group were not detected until 
genealogies were checked later. They are indicated in Table I and are 
included in the results because they show no peculiarities. 

Frazer (°33) has pointed out that ossification times differ from 
people to people and place to place, and are subject to the influence of 


*The material upon which this work is based was collected on a University of 
Adelaide anthropological expedition financed largely by grants from the Wenner- 
Gren Foundation for Anthropological Research, Incorporated, New York, and the 
University of Adelaide. 











266 A. A. ABBIE AND W. R. ADEY 


a number of environmental variables. Whether there are, in addition, 
really significant differences between different ethnic groups still remains 
to be settled, however. 

It is not pretended that the picture presented here is characteristic 
of the aborigines as a whole. Even given a common hereditary ossifi- 
cation pattern there remains ample scope for substantial differences 
between these people of the relatively arid interior and those of, say, 
the more luxuriant north. The present series itself may not be homo- 
geneous from the dietary point of view. The group came under govern- 
ment care only recently. Consequently, while the older members passed 
their early childhood on a purely native diet, the younger have enjoyed, 
in addition, access to European food supplements, together with some 
measure of pre- and postnatal care (see Campbell and Barrett, *53). 

Disease may affect bony maturation (Harris, °33). This possibility 
is not considered here, mainly because of ignorance of the local endemic 
diseases. Yaws (frambesia) must be fairly common, even in youth, 
judging from the incidence of “ boomerang” or “saber” tibia (plate 
1, A), but its effects on growth are unknown. There was no evidence 
in such people of interference with ossification times. 

The number examined were relatively small and age estimations were 
sometimes uncertain. Consequently, the findings are presented with 
some reserve. It is hoped that they will stimulate the prosecution of 
more extensive and better equipped investigations on as many other 
tribal groups as is still possible. 


MATERIAL AND METHODS 


The x-ray apparatus employed was a Watson Victor M 1 portable. 
Power was supplied by a petrol-driven field generator and a portable 
transformer. The films were 8 inch < 10 inch “ Kodarex ” which require 
no intensifying screen. Exposures were mostly at 30 kv. and 10 to 12 
milliamps for 2 to 5 seconds, at a distance of 30 inches. Apart from a 
field check, development of the films had to be deferred until return to 
Adelaide but the results were almost uniformly satisfactory. 

Because of limitations of time and material, radiographs were re- 
stricted to the wrist-hand and ankle-foot regions. Both sides were 
radiographed simultaneously. All the subjects also underwent a full 
anthropological and dental examination. 

Age estimations presented some difficulty. Only in the case of the 


very young were there any reasonable records. Even then there was 
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TABLE IL: AGE INCIDENCE 
MALES FEMALES 
Age Group | Ngalia Other Ngalia | Other 
< 1 mth - _ 1 (Sweeks) - 
1-6 mth. - _ 1 @mth.) - 
6-12 mth. | 1 (6mth) - = = 
12-18 mth. | 4 (i7mth) _ 1 (16mth) - 
18-24 mth. | 4 (20mth) - - - 
2-3 yr 1 (29mth) _ - _ 
3 yr. 3 - - _ 
3-4 yr. - - 1 - 
~ -& yt! - - { - 
5 yr | 1 _ 2 - 
5-6 yr.| - - 1 - 
6 yr { - 2 - 
6-7 yr; | 1 (Amanda y 1 - 
7 ye! 5 - 4 - 
7-8 yr. l 1 - - - 
& yr! 2 - 2 _ 
8-9 yr. _ -_ 2 _ 
oy) | t (eee) - 
10-41 yr. | - _ 1 1 (Pintubi) 
{1 yr. | 1 = — - 
41-12 yr. | 1 a _ _ 
{2 vr. 1 { ( — ) -_ _ 
1S yr. ; -_ 0 { - 
wn“ yt. 1 - 1 - 
14-415 yr. = { — 
15 yr. - - - - 
16 yr 1 _ 2 1 (Pintubi) 
17 yr. - -_ 1 - 
18 yr - - { - 
> 19 yr. - - { - 
Totals 24 3 29 2 




















* This lad was very precocious physically. 
Enclosed group is that used for study of Foot phalangeal epiphyses 
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room for an error of several days since, frequently, babies are born while 
their mothers are on bush excursions, and memory for dates soon fades, 
All age estimations were checked by such external evidence—mainly 
physical, and particularly dental—as could be secured. Even with the 
greatest care, however, there sometimes remained room for error of up 
to some months. 

The excellent kindergarten and school organization on the settlement 
secured a maximum attendance of the available children from 3 to 1] 
years of age. On the other hand, there is a regrettable shortage at both 
ends of the series. The youngest female was 3 weeks of age but the 
youngest male was 8 months. Obviously, it is misleading to regard 
these two ages as comparable, but when a center appears before 3 weeks 
they must be put in anachronistic juxtaposition when trying to establish 
a terminus a quo. There were sufficient adolescent females to provide 
a reasonable terminus ad quem for that sex but most of the adolescent 
males were out hunting or away in employment, leaving only a very 
precocious lad of 14 years and another—probably more typical physi- 
cally—of 16 years. 

In dating the results it must be remembered that if a center is 
already present it must have appeared at an earlier stage: this is indi- 
cated by the sign < before the age given. Only occasionally was a 
center caught at a rudimentary stage that could fairly be assumed as 
close to its first appearance. If any epiphysis had not fused at the 
maximum age in the series it was assumed that it would fuse later: 
this is indicated by the sign > before the age. In the case of epiphyses 
of metacarpals, metatarsals and phalanges it was possible to determine 
before complete fusion whether the secondary centers were proximal or 
distal. After fusion, of course, such distinction is sheer assumption as 
is, indeed, the postulation of any epiphysis at all. Complete fusion was 
recorded only when all trace of an epiphyseal line had disappeared. This 
is later than the standard adopted by Greulich and Pyle (’50) and may 
set back the apparent time of fusion. However, it seems the only way 
of being certain that fusion is complete. 

The series comprised 27 males, including 3 of mixed parentage, and 
31 females of whom 2 were pure Pintubi (a neighboring tribe). Their 


age incidence is shown in Table I. 
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, The male child of 20 months showed this epiphysis, but on the Left hand only. 
A = appearance. F = fusion. 


rile 
les, 
wae 4 SUPERIOR EXTREMITY 
the | TABLE I]: Distal Radial and Ulnar Epiphyses 
up Appearance Fusion Ulnar__Styloid 
Bone M F M E M E 

ent Radius X20 mths.|< 16 mths.| > 16 yrs. |18-19 yrs.) — — 
i | Uina KS yrs. |X S yes. | >l6 yrs. [18-19 yrs./6-7 yrs. |6-9 yrs. 
oth | 
7 TABLE Il: Carpus 
7 Lunate | Scaphoid] Trapezium| Trapezoid 
“ C- pitate | Hamate |Triquetrum| Semilunar| Navicular Mult. Maj. | Mult. Min. |Pisiform 
~ mie [Mle [Mie |[Mi|FI|MiF i M|F |M|F IME 
ride Tels | 8 k4 Ke0kie <5 (<4 <5 [<4 [<5 [<4 [<5 KS [e140 
ent mths} wks.| mths.jmths.| mths./mths.| yrs. | yrs. | yrs.| yrs. | ys.) yrs, | yes. |yts. | yrs.) yrs. 
—_ “hook” tubercle ridge 
| r at at at 
ysl- | 12 | 11 -i2ju-12} 12 | M4 

| yrs.| yrs. yrs.| yts.| yts.| yrs. 
r is 
di. | TABLE IV: Metacarpal Epiphyses 
Sa I pag il p 4 
| as M F M F M F M F M F 
the App. | 20mths|<16 mths.|<20 mths|{16mths. K24yrs. |{16mths. 2% yrs. {16 mths yrs. |<!6mths. 
ter: Fuse |>i6yrs. |<19 yrs. <i6 yrs. KAD yrs. 16 yrs. |<19 yrs. [ib yrs. Ki9 yrs. (ib yrs. {19 yrs. 
yses 
7 TABLE V: Hand Phalangeal Epiphyses 
or Dig. 1 Dig.2 Dig.3 Dig. 4 Dig.5 
‘= | MiF {[M{F |]MI{E Ml F M_[ F 
oes AS yes. <i6mihs|2omta Nomis 2Omths.{tGmths, 2mths mths. 2Omths mths 
alll | F Ki6yrs. Ki9yrs. [i6yes. Ki9yrs. Kb yrs. Ki9 yrs. KO yrs. [19 yrs. [6yrs. [A9yrs. 
all AK20mths. Kiomths.|<3 yes. KlOmths.|(20mths Ki6mths. (20mths.|6mths. |{3 yrs. |(l6mths. 
way 2 

F Ki4yrs. (i9yrs. Ki4yrs. M9 yrs. M6 yes. M9 yrs. Myers. |i9yrs. 4yrs. [19 yrs. 
and 34 7 _ [yrs Komths.K2'2yrs.Ki6mths. 2 eyes. |émiths. {S yrs. i6mths. 
heir F id yes. 19 yrs. M4 yrs. (9 yrs. 1414 yes. 9 yes. [44 yrs. |(19 yrs. 

| 
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RESULTS 

Superior Extremity 

Distal radial and ulnar epiphyses (Table IT) 
Radius In the male the distal radial epiphysis had not appeared at 
17 months (single case) but was present at 20 months (single case), 
Fusion was not complete in the lad of 16 years. In the female this 
epiphysis was bilaterally present at 16 months (single case). Fusion 
was incomplete at 18 years but complete in the woman over 19 years, 


Ulna In the male there was no trace of a distal ulnar epiphysis at 
3 years. The series did not include a 4 year old male but the epiphysis 
was well developed bilaterally at 5 years. In the female, the center 
was absent at 4 years, present at 5 years. Fusion appears to be about 
the same time as for the distal radial epiphysis. 

These times are not markedly different from those found in 
Europeans. 


Carpus (Table III) 

The order of ossification of these bones betrays no peculiarities 
although there is always some inter-variation between trapezium and 
trapezoid. In the present series one case (No. 2: Pintubi & Ngalia male, 
aged 9 years) showed apparent retardation on the radial side (right 
hand only), with the trapezoid ahead of the trapezium (plate 1,(C). 
Another case deserves mention. This (No. 56: a Ngalia male, aged 
3 years) showed bilaterally a peculiar center which could be interpreted 
as either a proximal, additional, epiphysis for the second metacarpal or 
as a slightly ectopic—and precocious—center for the trapezoid (plate 
1, B). 

In all instances the female appears to be ahead of the male, both 
in the appearance of a center and in its maturation as indicated by 
“ facetting,” and the development of special characters. Further, while 
in the male the pisiform center appeared at about the same time as 
the sesamoids opposite the head of the first metacarpal, in the female 
the pisiform was 3-4 years ahead of the sesamoids. 

No other peculiarities in the way of extra or ectopic centers or 
bipartite bones were noted. The times of appearance of most of these 
centers lie within the recorded European range. However, the tri- 
quetrum is earlier, as is the appearance of carpal maturation when 
compared with the figures given by Greulich and Pyle (50). 
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Metacarpal Epiphyses (Table IV) 

All primary centers were present in the youngest child examined, 
a female aged 3 weeks. Epiphyses for the first metacarpals were all 
proximal, those for the remainder all distal—with the possible exception 
of the second (No. 56) already discussed. There was no trace of any 
supernumerary digit or other variation and the only sesamoids seen 
were those opposite the head of the first metacarpal (flexor-abductor and 
adductor sesamoids) . 

The time of appearance of these centers in the female is consistently 
earlier than in the male. Against most European standards (Greulich 
and Pyle give generally earlier times) all the epiphyses in the female 
and those of the first and second metacarpals in the male are precocious. 
Male time of fusion here appears ahead of European, female about the 
same as European, but the instances are too few for assurance. 


Hand Phalangeal Epiphyses (Table V) 

All the primary centers were present in the youngest child examined. 
All the epiphyses observed were proximal; there was none extra or 
ectopic. In both sexes in both hands the epiphysis of the second 
phalanx of the thumb was earlier and always larger than that for the 
first phalanx (plate 1, B,C). 

The female adolescent series gave the general order of fusion as 
distal: proximal: middle. The males were too few to supply much 
information in this respect. As with the metacarpals, the female seems 
to be ahead in appearance and behind in fusion, but the series is 
restricted. 

The times of appearance of these epiphyses lie within the European 
range as given by Greulich and Pyle. The times of fusion may be earlier 
in the males and later in the females, as compared with Europeans. 


Inferior Extremity 
All radiographs of the foot-ankle region were taken in the antero- 
posterior position with the foot plantar-flexed. This secured the maxi- 
mum “coverage” per film, but in some instances the lower end of the 
tibia and fibula overlapped the hinder tarsal bones and obscured the 


picture in that region. 


Distal Tibial and Fibular Epiphyses (Table VI) 


Tibia In the male the lower epiphysis was already present in the 
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yé 
INFERIOR EXTREMITY bar 
TABLE NI: Distal Tibial and Fibular Epiphyses. “ 
Appearance Fusion | 
Bone M F M F ; 
Tibia <8 mths. <4 mths. |<16 yrs. |<16 yrs. ; 
Fibula |<20mths. 3-4 yrs. |<16 yrs. |<16 yrs. ‘a 
* May have appeared at 16 months, but that plate was obscured in this region. | th 
as 
TABLE VII: Tarsus. 
Calcaneum} Talus Cuboid | Cuneiform [Cuneiform | Cuneiform |Navicular bo 
Os Calcis |Astragalus na Scaphoid M 
Mit [MIF I|M/F [MIF IMIEIMIFEIMIE 
C813 KS KS KS F 3K 8 K4 K20 K3-41K29 K3-4K29 [3-4 I 
mths] wks./mths.|wks.i mths. lwks./mths.imths|mths.| yes. imths.| yes.imths.| yrs. 
Per. Long 
7 ate Tubercle al 
be 7 5 th 
po yrs. yrs.i yrs. In 
* Just appearing. fae 
wi 
TABLE Vill: Metatarsal Epiphyses. is 
I i juny WW y by 
M F M F M F M F M F 18 
App.j< Syrs.K3-4yrs.K 3 yrs. K3-4yrs.K3 yrs. K3-4yts./* Syrs.K3-4yrs.K Syes. |3-4yrs. dit 
FuseKl4 yrs. (19 yrs. |<16 yrs. [18 yrs. Ki4 yrs. K16 yes. |<14 yrs. KiB yrs. Ki4 yrs. KiB yrs. 
* Just appearing. | tal 
Note: Metatarsal V styloid. Male at6-8yrs. Female at 5-7 yrs. no 
TABLE IX: Foot Phalangeal Epiphyses. | ye 
Dig. 1 Dig.2 Di¢g.3 Dig.4 Di¢.5 M. 
| =| M F M F M F M F M_| F 
' A} 29mths] {6mth.K20 mth. Ki6mth. K20 mth.(16 mth. K20 mth. K16 mth.|<3 yrs. 6 mth. 
F Kt6yrs. Ki6 yrs. |i6 yrs. P19 yrs.| f6yrs Ki9 yrs. | 16yrs. X19 yrs.| 16 yrs. K19 yrs. | ™ 
oA 2Omth.K16 mth. K20 mth. 3-4 yrs.K20 mth.K4 yrs. KSyrs. |?(3-4yrs|? 3yrs.|?(3A4yrs | = 
F /16 yrs. [144]yrs| 14 yrs. |M61Byrs.| 14yrs. Ki6yrs. | 14yrs. 6 yrs. (i4 yes. K (6yrs. ob 
34 < i Syrs.|(3-4yrs.|(3 yrs. K5-4yrs./(Syrs.3-A4yrs|? 3yrs./2<s-dyrs pr 
F tAyrs.| 14-l6yts} 14yrs.Ki6 yrs.| fdyrs.|6yrs.| (14yrsi<i6yrs, 
A = Appearance. F = Fusion. “A 


ts. 








OSSIFICATION IN A CENTRAL AUSTRALIAN TRIBE 273 


youngest examined—8 months. In the female the center appeared 
between the third week and the fourth month. In both sexes time of 
fusion was about the same but was obscured by the overlapping referred 
to above. 


Fibula In the male the lower epiphysis was present before 20 months. 
In the female it may have been present at 16 months, but the radiograph 
for the individual of that age was obscured and the next—3 to 4 years— 
is the earliest that can be cited. Time of fusion appeared to be about 
the same for both sexes but estimates are subject to the same limitations 
as for the tibia. 

These figures show no clear-cut sexual distinction. In both sexes 
both bones fall within the maximum European range as given by 
Meschan (51 a, 


Tarsus (Table VIT) 

The fact that the major centers for the caleaneum and talus were 
already well present in the youngest child examined (3 weeks) supports 
the contention that these centers—as in Europeans—appear before birth. 
In the same child the condition of the cuboid was consistent with the 
view that its center appeared at or just after birth. The navicular falls 
within European extremes. The order of ossification of the cuneiforms 
is an accepted variant but the times are early. Maturation, as indicated 
by characteristic sculpturing and the appearance of distinctive features, 
is generally earlier in the male but otherwise there is little sexual 
differentiation. 

There was no evidence of a special center for an os trigonum on the 
talus, an epiphysis on the caleaneum or the tubercle of the navicular, 
nor were any unusual centers discovered. The female aged about 19 
years had bilateral tibialis posterior sesamoids. 


Metatarsal Epiphyses (Table VIIT) 

All primary centers were present in the youngest child examined 
(3 weeks). The epiphysis on the first metatarsal was always proximal 
that on the remainder always distal. No ectopic or extra centers were 
observed and there was no separate center for the tubercle (styloid 
process) of the fifth metatarsal. 

Taken all round, the female may be ahead of the male in times of 
appearance of these centers but there is little evidence except, perhaps, 
for the fifth metatarsal. The times are close to those given for Euro- 
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peans. Fusion appears to be earlier than in Europeans in the male, 
later in the female, but again the evidence is scanty. 


Foot Phalangeal Epiphyses (Table IX) 

All primary centers were present in the youngest child. So far as 
could be determined the epiphyses were proximal except for a possible 
variant in the little toe, mentioned below. In the big toe in both sexes— 
as in the thumb, but to a lesser extent—the epiphysis of the second 
phalanx appeared earlier, and was always larger than the epiphysis of 
the first phalanx (plate 1, E, F). 

Table IX gives times of appearance and fusion for all the epiphyses 
of all the toes but it should be emphasized that not a single foot attained 
this ideal. In both sexes in all feet both phalanges of the first toe, and 
the first phalanges of all the other toes, behaved in “ text-book ” fashion. 
However, the middle and distal phalanges of the second, third, fourth 
and fifth toes were subject to considerable variation. This is shown in 
Tables X and XI. Here the age group is chosen to ensure inclusion of 
all epiphyses likely to appear (Table I). The males (aged 5 to 12 years) 
and females (aged 4 to 13 years) are approximately equal in numbers. 

It will be noted that from the second toe laterally there is a pro- 
gressively greater tendency to suppression of epiphyses in the middle 
and distal phalanges (Table X). Either phalanx, or both, may be 
involved but the middle is affected more than twice as often as the 
distal. It is also apparent that the incidence is much greater in the 
females in whom more than half had all 4 lateral toes affected and 
more than three-quarters had at least three toes involved. The male, 
on the other hand, more frequently showed involvement of only one 
or two toes. 

In neither foot of either sex was there a little toe with its full 
complement of bony elements. Moreover, the bones present varied con- 
siderably in size and shape and often bore little resemblance to either 
primary or secondary centers (plate 1,D-L). Such bony fragments 
are best referred to as ossicles. The maximum possible number of such 
ossicles in the distal two joints would be 4: in fact, in no little toe 
were more than three observed, two was the commonest number, but 
one alone occurred occasionally (Table XI). There has evidently been 
some suppression (or, possibly, very premature fusion) of some centers, 
but exactly which is not always clear. Thus the middle of three ossicles 
may, apparently, behave like the epiphysis of the distal phalanx and 
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No. 53, Ngalia girl, aged 4 years 
Ni 3. Ngalia boy aged 9 years 
N 78, Ngalia girl, aged 16 years 


No. 42 
produced here 
rected during maturation.) 


Ngalia girl, aged 16 years 
istrate the pomt that deficiency of ossicles in 
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fuse with that or, more commonly, act as the diaphysis of the middle 


phalanx and treat the proximal ossicle as an epiphysis (plate 1, K,L), 
That, at least, is what examination of these radiographs suggests but a 
longitudinal series would be necessary to be sure in any one case, 

A point of interest is that these osseous changes in the feet were 


nearly always symmetrical. 


SUPPRESSION OF BONY ELEMENTS IN FEET 
of 18 Males, aged 5-12 yrs.and 19 Females, aged 4-I3yrs. 


TABLE X: Number of Toes Affected. 
5th Only [St & 4th | 5th gth| sth | 4th, aa 
& 3°4 | 34624) Uncertain | Totals 
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TABLE X1: Number of Phalanges Affected. 


























Suppression of Epiphyses 
2n4 Toe [54 Toe | 4 Toe |Ossicles} Sth Toe [Uncertain 
Phalanx} M/F [M|FIM/F fm |uMl] el MF 
a 2 0 5 
Distal 5 O 0 5 S 
Middle | © 3 9 2 12 2 
5 Is} | lof 1 2 
Both 0 7 5 1 ! 
0 1 6 1 












































* See text concerning 5™ toe. 


DISCUSSION 


Despite some unfortunate gaps in the evidence a few conclusions 
of interest emerge. In the first place, the times of ossification displayed 
by the aborigine lie within the total European range as given in 
standard works by Quain (715), Frazer (733), Buchanan (°49), Gray 
(49), Greulich and Pyle (50), Cunningham (7°51). Meschan (751), 


and others. The most that can be said for the aborigine is that he 
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usually falls into the earlier part of the range. This would be in line 
with Campbell’s (725) observations on the generally earlier dental 
maturation in aborigines. 

Most writers emphasize that in any ethnic group the female both 
begins and ends ossification earlier than the male. The figures for the 
aborigine generally support the contention for an earlier start in females, 
but what evidence there is suggests a later finish. Here, however, the 
facts are so scanty that it would be unwise to insist upon the reliability 
of this interpretation. 

One of the most striking results of this—admittedly limited—survey 
is the almost complete absence of a wide range of possible bony varia- 
tions as described for the European foot (see, e.g. Wood Jones, 744). 

In default of any positive evidence for specific ethnic peculiarities, 
attention may be directed to what seems to be the most interesting 
outcome of this study, namely, the evidence for “ degeneration ” or 
“recession ” on the lateral side of the foot. It is something of an article 
of faith that such a process is occurring in the European foot, but details 
are hard to find in standard works on anatomy and even in books devoted 
entirely to the foot. These all seem to contain an underlying assump- 
tion that the full number of bony elements is the rule and that any 
deficiencies of the moment are the result of delay, to be repaired in time. 
However, Flecker (°32) has shown that this is not necessarily the case. 
For urban children of mostly European parentage he reports a condition 
very similar to that described here for the aborigine, and with the same 
heavier incidence in the female. 

When it receives any attention at all, this suppression of osseous 
elements in Europeans is usually attributed to the effects of civilized 
existence and, particularly, footwear. The greater female incidence has 
been considered confirmatory evidence. But, the fact that the same 
phenomenon appears at an equal intensity, and with a similar sexual 
distribution, in aborigines who have never known civilized life or foot- 
wear exculpates both factors as possible agents. The process, then, must 
be part of a common human heritage, associated with increased emphasis 
upon the medial side of the foot in the interests of stability. The greater 
intensity in females may reflect special stresses imposed by the broader 
pelvis and associated bony adjustments. 


i 
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SUMMARY 


e This paper records the radiological condition of the wrist-hand, ankle- 
foot regions in 58 Central Australian aborigines of ages from 3 weeks 
to 19 years. 


~w 


The study as a whole disclosed no ethnic peculiarities. 

3. Ossification falls within the time range assigned to Europeans but 
usually towards the earlier rather than the later limit. 

4. An account is given of phalangeal “ degeneration ” in the foot similar 

to that found in Europeans. 
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FACTORS OF DOSAGE AND HOST DETERMINING 
ANTIBODY RESPONSE TO SECONDARY 
ANTIGEN STIMULUS * 


BY JOHANNES IPSEN, JR. 
Institute of Laboratories, Massachusetts Department of Public Health, 
SECONDARY antigen stimulus, according to Glenny and Siid- 

A mersen (°21) is defined as the injection of antigen into an actively 
immunized animal. The active immunization may be due to natural 
infection or to artificial immunization with the same specific antigen. 
The individual may or may not show evidence of such previous treat- 
ment by presence of specific antibody in the serum. The effect of a 
secondary stimulus is a distinctively higher immune response than that 
of a primary stimulus. Smaller doses of antigen are necessary to pro- 
duce a secondary response, than those needed for a primary response. 

The primary stimulus, however, is preferred in assays of the relative 
potency of toxoids, because greater uniformity of response is expected 
from a non-immune population. In field trials of toxoids in human 
population samples, one is often faced with a mixture of partly immune 
and non-immune individuals, so that the interpretation of the response 
is equivocal. Accurate information of previous immunization is often 
essential for a reduction of the variability through covariance. The 
following presentation is an attempt to evaluate the importance of 
measurable independent variables and their contribution to the variation 
in secondary response in human beings. 

Response is in this work measured by the antibody titer which is one 


measurement of immune response. Other response measurements such 


‘The experiments were conducted under the sponsorship of the Commission 
on Immunization, Armed Forces Epidemiological Board, and were supported in 
part under contract with the Office of the Surgeon General, Department of the 
Army. Presented at the Third International Conference on Biometrics, Bellagio, 


Italy, September 2, 1953. 
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as percentage of survivals or percentage of skin reactions to antigen or 
toxin, are less informative since they have the limitations of all-or-non- 
reactions. A fundamental relationship between antitoxin titer and all- 
or-non-reactions has been well substantiated recently by Holt (753), 
and Cavalli-Sforza (’47). 

The discussion is limited to the antitoxin titer at a certain time after 
the antigen injection, which approximates the time of maximum response, 
The quality of the antitoxin which recently was subjected to intense study 
by Jerne (°51) is also omitted from this discussion, where antitoxin titer 
simply denotes the highest dilution of a serum which neutralizes a given 
dose of toxin to a certain end-point reaction. Since the titers are 
uniformly determined throughout the assay, this study deals with com- 
parison of response, where the absolute value of the response is unimpor- 
tant. The notation of the titers in Antitoxin Units per ml (AU) is in 
conformity with present convention. 

Titrations were performed in mice with constant toxin doses and 
serial 1:2 dilutions of serum. 

For a complete analysis of the factors determining a given secondary 
response the following variables must be included in the experiment and 


each must be varied over a sufficiently wide range. 


1. (y) The dependent variable which is the titer found after the 
secondary stimulus. Use of logarithm increases the homosche- 
dasticity of the system somewhat, but there is a tendency of larger 
error at smaller antitoxin levels. 
2. The independent variables are: 
a. (z) The size and potency of the secondary stimulus, also 
transformed into logarithms. 
b. (z.) Intensity of previous immunization, which includes 
the dose and potency of earlier antigen exposure and the time 
elapsed since it was administered. 
ce. (v) Pre-immunization titer, i. e., the antitoxin titer of the 
individual immediately before the secondary dose is given. 
d. (z,) The inherent ability of the individual to react to an 
antigen stimulus. This factor can only be assigned quali- 
tatively to groups of subjects with given species characteristics 
or some somatic abnormality. 


The immunization history, the pre-immunization titer and the in- 
herent immunizability together form the immune status. Rarely all the 
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information is available in a human field study. Accurate knowledge of 
the time and size of previous injection is often missing, so that sub- 
divisions such as “ number of injections,” “time since last injection,” 
ete., are the only measurable qualities. The pre-immunization titer is 
sometimes considered a unique measure of the immunization history. 
Although a high titer suggests effective previous immunization, low titers 
or titers below the measurable limits may represent all sorts of immune 
status including the so-called “ non-immune ” where an antigen injection 
would provoke a primary response only. The uncertainty as to whether 
a person has primary immunity is larger in case of diphtheria than of 
tetanus, since natural immunity is possible in diphtheria while artificial 
immunization is the only possible way of acquiring tetanus immunity. 

Immunizability is measurable as a group phenomenon rather than as 
an individual characteristic. It can be shown that animal strains of 
various genetic composition require different doses of antigen to provoke 
primary response. Thus, we have been able to show that two pure inbred 
strains of mice differ 12 times in the amount of tetanus toxoid necessary 
to induce 50% primary immunity to challenge with tetanus toxin. 
Similar differences were noted in guinea-pigs with respect to diphtheria 
immunizability (Prigge, °35; Scheibel, °43). 

The influence of inherent immunizability on secondary response needs 
systematic research. The experimental data to be presented afford an 
opportunity to analyze the four components of primary dose, pre-immu- 
nization titer, secondary dose and inherent immunizability in a group 
of young men. It was possible to conclude that primary dose and primary 
response were interchangeable measurements and that immunizability is 
a property which is more pronounced in secondary response than in 
primary response. 


EXPERIMENTAL DATA 


The experiment was originally designed to assay the relative potency 
of four aluminum precipitated toxoids. The subjects were inmates of 
the Wrentham State School, Massachusetts, of which 128 young men, 
aged 15-25, were selected because of no previous history of tetanus 
immunization, which was furthermore confirmed by the fact that none 
of tie subjects showed measurable antitoxin in the blood before the test. 
Each of the four toxoids were given in either of two doses, 0.1 ml and 
0.5 ml, respectively. These eight treatments were repeated 4 weeks after 
the first injection, but the individuals were distributed according to an 
8 X 8 latin square with the primary injections forming one axis, and the 
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secondary injection forming the other axis. Thus all 64 combinations 
of primary and secondary injections were given. (See table 1.) 

The columns of table 1 contain the antitoxin titer (in logarithms 
plus 3). In the left hand part of each column (v) are those titers deter- 
mined 28 days after the first injection; in the right hand part are the 


titers 14 days after the secondary injection (y). 


“ee 8 £2400 « 


1125 AU 


De 


8s less than lower limit of titration (0.0¢ 


0.1 indicates value 








RELATIVE POTENCY OF DOSAGE 


A previous report on these data, (Ipsen, *53A), dealt with the esti- 
mate of the relative potency of these four toxoids, which had also been 


assayed in both mice and guinea-pigs. The potency estimates were 


obtained from this assay in humans by construction of reaction curves 


relating primary response (v) to dose of toxoid. 
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SACTION CURVES OF Four TETANUS TOXOIDS IN HUMANS 
(PRIMARY RESPONSE) 


Fic. 1. RE 


The reaction curves clearly place Toxoid B as the most potent, 
Toxoid C of average potency, and A and E as the least potent. The 
same potency relationship was also found to hold when the secondary 
responses were related to the secondary dose. 

In this particular design each of the groups which received a secon- 
dary dose consisted of individuals who had the same average primary 
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immunity status. Hence, no correction for primary immunity is neces. 
sary for the estimation of relative potency based on secondary response, 

The eight treatments can now be placed in order of relative potency 
as indicated in table 2. It should be noted that the potency order of 





























the same four toxoids was found to be identical to that determined jp 
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TABLE 2 ce, 
Mean Primary and Secondary Response, by Potency of Toxoid Dose ' M 
Jif 
Relative Log Rel. Primary Secondary Primary | 
Toxoid Dose Potency Potency Number Response Number’ Response _ Response’ | 
E 0.1 ml 1.0 0.0 16 0.100 16 1.79 0.51 — 
A 0.1 ml 1.6 0.2 16 0.119 15 2.27 0.40 & 
C 0.1 ml 3.2 0.5 16 0.213 14 2.49 0.66 | v 
E 0.5 ml 5.0 0.7 16 0.400 14 2.16 0.36 2v 
A 0.5 ml 8.0 0.9 16 0.475 18 2.48 0.43 ‘ 
B 0.1 ml 8.0 0.9 16 0.550 16 2.61 0.57 : 
C 0.5 ml 16.0 1.2 16 0.888 17 2.84 0.51 ess 
B 0.5 ml 40.0 1.6 16 1.150 18 3.25 0.47 y 
RTA Stee ib decals Ul 1 (1) ans a re Mi 


* Mean titer of group before secondary dose was injected. 
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REGRESSION ANALYSIS OF PRIMARY AND SECONDARY DOSAGE, 
DISREGARDING IMMUNIZABILITY 


For the biometric analysis, we shall use y (log. secondary antitoxin 
titer) as the dependent variable, and let the above logarithmic relative 
potencies substitute for x, and z, as measurements of the primary and 
secondary dosage. 

Before assumptions can be made as to primary immunizability, a 
multiple regression analysis was performed with linear and quadratic 


terms of x, and z and their interactions. Only three terms were found 


TABLE 3 


Matrix of Sums of Squares and 
Sums of Products of Deviations from Mean 


128 Observations 


z Bs ZL, @, 22, y 

“s 31.56 7 ieee 94.91 
@s 0.12 30.40 11.25 
22, 42.39 23.04 42.39 22.24 
2, 2.95 1.80 1.82 24.00 12.31 
2, 8.12 0.38 5.17 17.68 20.39 —2.24 
y a 100.2437 
Means 0.778 0.750 0.580 0.125 —0.120 2.507 

2 v 2v a, 22, y 
2 31.56 win = 24.91 = 
1 0.21 46.11 12.79 
20 15.07 34.21 42.42 15.02 
@, 2.95 7.39 4.13 24.00 12.31 
za, 8.12 5.60 2.20 17.68 20.39 —2.24 
Means 0.778 0.487 0.377 0.125 0.120 2.507. 
50 Observations (v larger than 0.1) 
2 v 2v av, 22, y 

2 12.83 2 5.03 
v 1.59 16.25 0.96 
20 12.07 10.96 24.13 3.88 
a, 1.59 3.38 1.01 7.92 1.69 
2a, 1.65 3.81 0.86 7.45 8.39 0.40 
y 18.3442 


Means 0.806 1.090 0.847 —0.040 —0.064 2.746 
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TABLE 4 


Analysis of Variance 








Primary Dose (a,), Secondary Dose (z) and Immunizability (2,) on ! 
dev! 
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significant (see table 4, Analysis of Variance) so that a preliminary 
expectancy equation was computed. 


y = 1.243 + 0.8702, + 1.2862 — 0.667422 - « (1) 


on the basis of the matrix of sums of squares and sums of products of 
deviations. (see table 3). 

The fit of this expectancy for the overall sample of 128 individuals 
was next tested for smaller groups. Since the inmates of the institution 
represented a great variety of constitutions, it was pertinent to investigate 
whether grouping of the individuals in categories of their somatic and 
mental diagnosis would reveal any group which deviated in response 
from the overall expectancy to suggest abnormal immunizability of the 
group. 


TABLE 5 
Observed and Expected Means for Ten Diagnostic Categories. 
Expectation from Equation (1), Disregarding 
Inherent Immunizability, and from Equation (2), Including Group 


Immunizability 


Equation (1) Equation (2) 

DIAGNOSTIC Number Observed Expected * Student’s ” Expected es Student’s ” 

CATEGORY n Mean Mean e” Mean ae 
Undiff. imbecility 44 2.49 2.50 —(.04 2.54 —0.43 
Familial morons 31 2.65 2.55 +0.70 2.59 +0.43 
Central nerv. syst 

malformations 5 2.90 2.76 +0.41 2.77 +0.41 
Epilepsy and other 

defects 6 2.60 2.27 +1.03 2.36 +0.83 
Psychosis and mental 

deficiency 6 2.20 2.52 1.03 2.56 —1.26 
Post-traum. idiocy 6 2.82 2.66 +-0.49 2.72 +0.32 
Post-infect. idiocy 4 2.38 2.49 —0.30 2.55 —0.49 

102 2.56 2.52 +0.51 2.56 —0.09 
Mongoloids 21 2.08 2.49 2.04 2.07 + 0.05 
Dysendocrines 5 3.18 2.23 +2.52 3.15 +0.11 
* 42 (Obs Exp 2n e720 — (Obs Exp) *n 


0.5907 0.5011 
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Table 5 shows the observed and expected group mean for ten major 
categories which were classified by the medical staff of the institution, 
before the immunization program started. 
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Fie. 3. OBSERVED RESPONSES OF 21 MONGOLOIDS AND 
5 “ DYSENDOCRINES ” IN RELATION TO VALUES 
EXPECTED FROM EQUATION | 


While the observed titers of eight categories, composing 102 persons, 
agree well with expectancy equation (1), there are two groups which 
deviated to the borderline of significance. Mongoloids (21 persons) 
have a too low response, and “ dysendocrines” (5 persons) show an 
excessive response. 

A further analysis of the deviations from expected values within the 
two groups show that deviations are largest from low expectancies, while 


high expectancies give acceptable fits. 
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That mongoloids have less immune response than other individuals 
is evidenced in studies by Siegel (’48). The group “ dysendocrines ” is 
not a well defined one but since it seems to represent an immunological 
antipode to the mongoloids, it is worth including it in an analysis of 
immunizability. 

Since the two abnormal groups have approximately the same numer- 
ical deviation, the immunizability can be tested as a linear regression 


by assigning the factorial 7, = -+ 1 to the dysendocrines, 7, —=—1 to 
the mongoloids, and x, —0 to the remainder categories. The matrix 


of table 3 contains the significant components; others were also tested, 
such as the interaction 2,2, without giving additional information. 
The significance of the factor z, and its interaction with z is brought 
out in the analysis of variance. Table 4. 

With the solutions of the regression coefficients the expectancy equa- 
tion now takes the form: 


y = 1.549 + 0.6362, + 1.0012 — 0.435722 + 1.1977, — 0.872a,2--- (2) 
Equation (2) can be rearranged in the form: 
y = 1.549 + 0.636(1.8827, + zr.) + 1.0012 — 0.4352(2.00z, + 22): - - (3) 


which indicates that the inherent factor z, has about equal additive 
effects to the dual functions of the primary dose z,: one of increasing 
the secondary response, the other of decreasing the regression of the 
secondary dose z on the response. 

The different contribution of the primary dose and of the immuniz- 
ability is clearly shown in Figure 4 which presents the observed data in 
relation to the expected regression lines of y and z for four groups. 


The upper line represents the expectancy for the 5 dysendocrines who 


received an average primary dose of z, — 0.76, which inserted in equation 
(2) together with z, —1, gives 
y = 3.26 — 0.2002: - - (4) 
The lowest line represents the equation for the mongoloids, where 
zt, = 0.67 and 7, = 1 leading to 
y = 0.719 + 2.1622: - - (5) 


The two intermediate lines are the expectancy for two primary dose 
levels of the 102 “normals,” of which the upper is the expectancy for 
47 individuals with relative high primary dosage (zr. — 1.14) and the 
lowest line represents 55 persons with a lower primary dose (rz, — 0.33). 
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The equations are, respectively : 


y = 2.274 + 0.5052: - - (6) 
and 
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Fie. 4. SrcoNDARY RESPONSE OF THE VARIOUS IMMUNITY LEVELS, 
RELATED TO SECONDARY DOSE 


The difference between the two “normal” groups is much less 
pronounced than the difference between the mongoloids and dysen- 
docrines who had approximately the same primary dosage. The points 
represent the mean of 1-8 individuals of the same kind and with same 
secondary dosage. 
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DIFFERENCE IN IMMUNIZABILITY AT PRIMARY STIMULUS 


The difference in the two somatic groups amount to about 2 < 1.88 
— 3.76 in terms of log. primary dose as estimated from the secondary 
response (equation (3)). This means that the mongoloids would need 
about 5600 times more potent primary stimulus than the dysendocrines 
in order to obtain the same responsiveness to secondary stimulus. 

This difference is in contrast to the difference in inherent immuniz- 
ability which can be estimated from the primary response. 

The following data are available for the response (v) of the three 
groups after the first injection, and are presented in table 6. 


TABLE 6 


Primary Response in Three Groups of Different Immunizability 


Average Percent with Mean Primary 
Group Number Primary Dosage Measurable Antitoxin Response (0) 
“ Normals ” 102 0.77 41.2% 0.521 
Mongoloids 21 0.67 23.8% 0.229 
Dy sendocrines 5 0.76 60.0% 0.880 
Difference in Immunizability 
Mongoloids vs. Dysendocrines 0.35 + 0.31 0.79 + 0.40 


(in terms of primary dosage) 


The difference 
cant, and even with the large error involved, the difference is nowhere 


n primary response is noticeable, but hardly signifi- 


near the extreme difference of 3.76 found from the secondary response of 
the two groups. 

This discrepancy between estimates of immunizability obtained by 
primary and by secondary response needs further study, which we are 
now pursuing in inbred mouse strains. (Ipsen, *53B). 


RELATION OF PRIMARY RESPONSE TO SECONDARY RESPONSE 


The foregoing analysis showed that information of both inherent 
immunizability and of the primary immunization is necessary to form 
estimates of the response to a certain secondary dose. 

The question now arises whether the primary response gives enough 
information of both the primary dose and the inherent immunizability, 
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so that the primary response can be used as an efficient covariance to 
the secondary response. 

The procedure used in the analysis was the same as used above, only 
the primary response (v) is used instead of the primary dosage (z,). 

Those primary titers which were below the limit of measurement 
were assigned the value 0.1 (as indicated in table 1). 

The matrix of sums necessary for the multiple regression analysis 
is found in table 3. The solutions of the regression coefficients and the 
means lead to the following expectancy equation: 


y = 1.747 + 0.567(v + 2.07z,) + 0.9172 — 0.4972(v + 1.742,)--- (8) 


This equation has almost the same form as equation (3) with insig- 
nificant difference between the coefficient of z, and v. It is particularly 
noteworthy that the coefficients of xr, are of the same magnitude in both 
equations. The analysis of variance (table 4, lower part) shows that z, 
and v give almost equivalent information, while the components of z, 
are of same order of magnitude. 

In other words, a given primary antitoxin titer reflects the amount 
of previous exposure to antigen but it does not predict the full extent of 
the response to a secondary stimulus. 

Since only 50 of 128 individuals had measurable antitoxin after the 
first dose, it is possible that the relationship expressed in equation (8) 
is biased because of the large number of v values which have been 
assessed the limit value of 0.1. 

In those persons where the primary titer can be measured, the effect 
of the inherent immunizability might be included in the primary titer. 
However, the equation obtained with these 50 individuals has the form: 


y = 1.61 + 0.72(v + 1.802,) + 1.272 — 0.822(v + 1.342,)--- (9) 


which is not significantly different from equation (7) obtained from all 
128 samples. 

The analysis of variance (table 4) shows that the components due 
to z, are at least as important as those due to v. Hence, it is concluded 
that inherent immunizability is a function which is independent of the 
primary response. 

DISCUSSION 
In the attempt to analyse the variables which influence the secondary 


response, the procedure has been restricted to analysis of linear com- 
ponents. The true functional equations are probably exponential func- 
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tions with curves tending toward asymptotes at high values of the 
respective variables. Such functions are not amenable to explicit regres- 
sion analysis, and no increase in information was obtained by adding 
higher terms of a polynomial expansion. The three measurable variables: 


(z) potency of secondary dosage 
(x.) potency of primary dosage 
(v) antitoxin before secondary immunization 


were all found to have some relation to the secondary response. The 
correlation, however, vanishes at high values of all or either of the 
variables. 

[t was further possible to recognize two somatic groups of individuals 
which deviated from the general expectancy in positive and negative 
directions, respectively. This discrepancy is ascribed to inherent factors, 
termed immunizability. While the primary response barely disclosed 
deviating immunizability, the factor was significantly noticeable after 
secondary stimulus. 

Because of the inhomogeneous character of the population used in 
this study it is highly likely that other individuals than the dysendocrines 
and the mongoloids have deviating inherent immunizability; however, 
these groups were the only ones which could be separated a priori on 
the basis of the available information. Other categorization, such as 
blood type, was tried with no result. 

Bio-assays based on secondary response are feasible, but they meet 


more obstacles than assays based on primary response. 


The restrictions are: 
1. Small secondary doses only will give significant regression of 
response on dose. What should be classed “small dose ” depends on 
the animal species used. In men, such doses are of reasonable mag- 
nitude, which appears from this study of tetanus toxoid and from 
that of Edsall, Banton and Wheeler (°51) who worked with small 


secondary doses of diphtheria toxoid. 


2. Individuals with high immunity status give less significant re- 


sponse, i. e., rise in titer after the secondary dose. 


3. Individuals with no measurable immunity may be of two cate- 
gories, previously immunized or not immunized. These inseparable 
categories yield widely different results, and immunization histories 
are necessary for efficient interpretation. 
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4. Inhomogeneity of the group with regard to immunizability intro- 
duces an important variable, which can rarely be determined a priori, 
except in pure inbred animal strains. This variable is most distinct 
after secondary stimulus. 
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INBREEDING COEFFICIENTS OF THE RAMAH 
NAVAHO POPULATION 
BY J. N. SPUHLER AND CLYDE KLUCKHOHN 
Institute of Human Biology and Department of Anthropology, 


University of Michigan; Departments of Anthropology and of Social 


Relations, Harvard University 
INTRODUCTION 


NY discussion of inbreeding should distinguish between two senses 
in which the term has been used: the social, and the biological. 


s, to mating 


Socially speaking, we may refer to endogamic mating, that 
which occurs within the limits of a socially-defined category, such as a 
community, a clan, a caste, or a class. Biologically speaking, inbreeding 
refers to the mating of individuals with one or more common biological 
ancestors. Inbred individuals are offspring of genetic relatives and 
inbreeding is the mating of genetically related individuals. 

The genetic consequences of inbreeding and endogamy may differ. 
Relatively small, endogamous and isolated human populations are not 
necessarily highly inbred. The degree of inbreeding in many natural 
and enduring local populations is slight compared to that characteristic 
of matings typical for some domestic and experimental animals. 

Throughout the remainder of this paper we shall be concerned with 
inbreeding in the biological sense, and not, except as noted, with endog- 
amy in the social sense. 

This paper describes the types, and gives a measure of the degree, 
of inbreeding in the Ramah Navaho population as of September first 
1948. A number of estimates of inbreeding coefficients based on various 
sorts of cousin marriages are available for human populations (see the 
summary in Neel, et al., 49). Measurements of inbreeding for relatives 
of all degrees, for a total human population seemingly have not been 
published heretofore. 
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INBREEDING COEFFICIENTS 


Wright (’22) developed a coefficient of inbreeding that gives (for a 
population) the departure toward complete homozygosity from the hetero- 
zygosity expected under a system of random mating. This coefficient is 
a direct measure of the average effects on zygosity due to inbreeding 
expected from the pattern of mating in any given pedigree, population, 
or population segment. Wright’s coefficient has the great advantage (not 
present in others, e.g., those of Pearl, °13, nor Krizenicky, *23) that 
the degree of inbreeding in an observed genealogy, however complex or 
irregular the pattern of mating, may be compared to that observed or 
expected under any regular system of inbreeding. In this way theoretical 
models, and the results of experiments, become useful for interpretation 
of human population data. 

In analysis of the Ramah Navaho data, we found it convenient to use 
the following definition of the inbreeding coefficient, i, identical with 
that of Wright (i,=f,): 


to = &(1/2)"7°(1— ig), 


where n is the number of steps (the path from parent to offspring) in 
a loop of descent closing on an inbred individual 0, Ls is the 1 of a 
common ancestor A, and summation is over all loops terminating in 0. 
If A is not inbred, the coefficient is reduced to 


ip = 3 (1/2)"". 


The quantity i, may be considered as the probability that an inbred 
individual with one given autosomal gene will possess two genes of 
common origin at the given locus. Also, 1) may be defined as the most 
likely fraction of all gene pairs to consist of homozygous pairs. It is 
important to emphasize that the coefficient as used here is relative to 
the average zygosity of a random-bred population. If p is the frequency 
of a gene, and (1—p) —4q is the frequency of its allele, the average 
heterozygosity of a population for the pair of alleles at equilibrium will 
be 2pq. The value of p and gq, of course, may vary with different gene 
pairs in the same population. The population effect of inbreeding on 
zygosity is to make the proportion of heterozygotes < 2pq, and thus, 
the proportion of homozygotes > p* + q’. 

Something of the genetic meaning of i is illustrated in figure 1 
where AA’ represents the genotype of a common ancestor. The prob- 
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ability that this ancestor will pass gene A to his child is 1/2. Likewise, 
the probability that he will pass gene A’ is 1/2. Now gene A may pass 
to 0 by the left side of the loop with a probability (1/2)*%, or by the 
right side with a probability (1/2)*. Thus the probability that 0 will 
be of homozygous genotype AA is (1/2)°. But gene A’ may pass over 
the same two paths (and 0 be of genotype A’A’) with equal probability— 
thus the factor 2 in the coefficient: 


to = 2D [(1/2)* + (1/2)?] — 2(1/2)°. 
Or, in more general terms, 


l 28° 
~ hm 





(1/2)™ + (1/2)*] = 23)(1/2)* = $(1/2)™"3 


where nz, is the number of steps in the left side of the inbreeding loop, 
nz the number in the right side, and n — nz + npg. 
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Fig. 1. DIAGRAMS TO ILLUSTRATE THE GENETIC BASIS OF THE INBREEDING 
COEFFICIENT. IN THE PEDIGREE ON THE LEFT INDIVIDUAL 0 Is INBRED: 
IN THE PEDIGREE ON THE RIGHT INDIVIDUAL 0’ Is Not INBRED 


The derivation and genetic significance of 1)—f is discussed in 
detail by Wright (°22, 51), Li (48), and Ludwig (°44). Wright derived 
the coefficient of inbreeding (f) by the method of path coefficients. 
Bernstein’s (°30) coefficient, a, is identical with f of Wright. Haldane 
and Moshinsky (39) give an algebraical derivation. Malécot (’48) has 
shown that the general formula for f may be obtained directly from 
probability theory. Using matrix methods, Fisher (’49) has proposed 
an inbreeding statistic which is identical with [— log,(1—f) ]. 

In the method of Haldane and Moshinsky, the relationship of fathers 
and mothers (parents of the inbred individuals in question) is classified 
for each inbreeding loop connecting the parents and the inbred individual 
0, and the value of % is found by such expressions as: for full r-th 
cousins s times removed 1) = 2-*"-**, and for half r-th cousins s times 
removed 1, 2-*r-s-8 the values of the i) being summed over all in- 


breeding loops that contain 0. 
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Calculation of the inbreeding coefficient is illustrated in figure 2 
which is the pedigree of the most highly inbred sibship of the Ramah 
Navaho population. In this pedigree 6 different inbreeding loops con- 
nect 0 and his 6 different common ancestors. The value of i, for each 
of the 6 loops is shown in the figure, their sum totaling .0977. This 





os 2 


Fie. 2. DIAGRAMS TO ILLUSTRATE THE CALCULATION OF THE INBREEDING 
COEFFICIENT. THE PEDIGREE OF THE Most INBRED SIBSHIP OF THE 
RAMAH NAVAHO POPULATION IS SHOWN ON THE Upper LEFT. THE 
INBRED INDIVIDUAL AND His 6 CoMMON ANCESTORS ARE SHOWN BY 
SoLtip Dots. THE 6 INBREEDING Loops CONNECTING THE INBRED OFF- 
SPRING AND HIs COMMON ANCESTORS ARE SHOWN BY HEAVY LINES IN 
THE REMAINING 6 DIAGRAMS 


value is a measure of the degree of inbreeding of individual 0. The 
value may be checked by classifying the relationship of 0’s father to 0’s 
mother and summing the values of i, for each type of relationship, 
following the method of Haldane and Moshinsky. The parents of 0 are 
related in three different ways: 1) full first cousins, i) —2-*, 2) full 
first cousins once removed, i, = 2-°, and 3) full third cousins, i, = 2°, 
the sum of these values being .0977. 





+ 26 + 
1G = 24 
1GIR = 2° 
36 = 2 

/ = 0.0977 
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The above remarks on the inbreeding coefficient assume that the 
alleles concerned are autosomal. Roughly 95% of human gene loci are 
autosomal, the remaining 5% being sex-linked, holandric, and partially 
sex-linked, (Spuhler, 48). Haldane and Moshinsky (’39) and Wright 
(51) have extended the definition of the inbreeding coefficient to cover 
these non-autosomal loci. Calculations of the inbreeding coefficient for 
sex-linked genes have been made for the Ramah population and will be 
presented in a later paper. The data of the present paper pertain to 


autosomal genes only. 


THE RAMAH NAVAHO POPULATION 


The Ramah Navaho are a local group of Navaho Indians residing 
in the vicinity of Ramah, west central New Mexico. The Navaho are 
an Athabascan-speaking people who have lived in the American South- 
west since at least 1,000 A.D. Unlike the Anglo-American, Spanish 
American and Pueblo peoples of the area, the Navaho do not concen- 
trate their population into villages. Rather, the community is composed 
of biological families (or more usually such families combined into 
extended families with matrilineal descent and matrilocal residence) 
widely spaced in groups of small houses or hogans. The Navaho sub- 
sistence pattern is a combination of sheep-husbandry and dry-farming 
agriculture. The closest Anglo-American towns are Gallup, 42 miles 
distant and about 11,000 in population, and Grants, 60 miles away and 
3,000 in population. 

Kluckhohn and Leighton (°46) give an introduction to the anthro- 
pology of the Navaho, with some emphasis on the Ramah people; 
Kluckhohn and Spencer (°40) provide a comprehensive bibliography of 
the Navaho literature up to 1940. 

The group at Ramah was founded about 1820 when a few families 
moved into the area to hunt and farm. They were preceded there by 
the village-dwelling Zuni Indians, but they occupied the area before the 
present Spanish American, Mormon and “ Texan” settlements were 
established. The Ramah group was reestablished soon after the Navaho 
returned from U. S. Army captivity at Fort Sumner, New Mexico, that 
is, in the decade 1870-1880. 

The Ramah population is isolated geographically from the larger 
concentrations of Navaho in northwestern New Mexico and northeastern 
Arizona. The Ramah group lives in an area of about 18 by 29 miles 
including some 505 square miles of which about 230 are Navaho-con- 
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trolled. Three other Navaho communities are in the general region, 
one 25 miles to the north, one about 50 miles to the east, and one about 
75 miles to the southeast. There has been some biological interchange 
between the Navaho local groups. 

On September first 1948 the Ramah group totaled 614 individuals, 
giving a Navaho population density of around 1.22 persons per square 
mile. Roughly 600 non-Navaho individuals occupy the same area 
{mostly Anglo-Americans, Spanish Americans, and Zuni Indians) with 
a total population density of about 2.40 persons per square mile (see, 
for brief additional information, Kluckhohn and Griffith, 751). 


CULTURAL FACTORS IN NAVAHO MATE SELECTION 


The explicit cultural regulations on Navaho marriage are mainly 
defined negatively. Thus in the case of the Ramah Navaho, one may 
not marry: 

1) a member of one’s own clan, 


2) a member of one’s father’s clan, 


-3 


3) a member of a clan “linked” to one’s own clan. 


a clan “linked ” to one’s father’s clan. 


4) a member 


The Ramah genealogies show 4 violations (out of 399 known mar- 
riages) of the first prohibition, 5 of the second. Transgression of the 
third and fourth regulations are more frequent but cannot be enumerated 
exactly because of disagreement among the Navaho themselves as to 
which of the clans should be considered “ linked.” In any case violations 
of these types are not as such biologically significant except perhaps at 
quite remote generations for which precise data are not available. 

The existence of matrilineal exogamous clans among the Ramah 
Navaho does make for outbreeding, both because marriage with first 
cousins on either side is specifically forbidden (but 4 sibships in the 
population have been produced by first cousin matings) and because, 
at Ramah, all members of the same clan (with a few exceptions) can be 
shown to be genetic relatives. In general, it may be said that the effect 
of clan and clan-group exogamy is toward a wide selection of mates 
within the local population and toward the procurement of mates from 
outside communities. It tends toward forcing genes from other Navaho 
groups, other Indian tribes, Spanish Americans, and Anglo-Americans 
to be spread widely (and, in the Ramah case, thinly) in the population. 
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(Of the 31% fertile matings shown in the genealogies, 35 or 11.1% 
involve a mate from outside the Ramah population.) To some extent 
also, the cultural interferences with mating tend to cancel one another. 

Clan membership may turn out to have other consequences for the 
genetic composition of this population which we are not prepared to 
demonstrate at present but which we hope to treat in later papers. We 
now call attention merely to the varying numerical composition of 
Ramah clans and to the fact that marriage between clans does not 
occur at random. 

Of individuals in the genealogies whose clan membership (or non- 
Navaho affiliation) is definitely ascertained, 77% are included within 
the 4 largest clans in this area, 97% are included within the 6 largest 
(see Kluckhohn and Griffith, 51, p. 405 ff., for a survey of the com- 
position as of September first 1948; data presented in that paper and 
in the present section overlap but also supplement one another). That 
is, all marriages of individuals from 15 additional clans plus those of a 
Yaqui Indian, a Walapai, a Laguna, a Spanish American, an Anglo- 
American, and 8 Chiricahua Apaches account for only 3% of the total 
marriages. 

TABLE 1 


3822 Ramah Navaho marriages classified by clan membership 


of the spouses 


Clan B C D -E F G H I. Total 
A 31 39 26 21 14 6 2 9 148 
B 95 31 10 8 2 1 6 83 
( 15 4 7 5 ~_ 16 47 
D 4 10 9 we 7 32 
E 8 l _ 3 12 
Total 31 64 7 4 #4 #216 #3 #41 «2322 ~~ 





9 


Of 399 known marriages (317 of these produced sibships and are 
analyzed for inbreeding in this paper) almost 88% are included within 
marriages of members of the 5 most numerous clans with each other and 
with members of 4 additional clans. Table 1 shows certain preferential 
tendencies for clans to exchange members in ways not determined by the 
system of linked clan exogamy. For example, 24% of the marriages of 
clan A members are with clan C members, whereas 14% is expected on 
a random basis in proportion to numbers of available adults. Preference 
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is even more striking in the case of the 16 matings of clan C with clan J 
members. The clan I column has another interest in that it is not a 
“native” Ramah clan. These marriages, except for two in the founding 
generation, represent the securing of mates from other Navaho groups, 

Additional cultural factors in mating, though not explicit regulations, 
are: permissible stepdaughter marriage (5 known cases in the history 
of this population) ; permissible polygyny with strong preference for 
sororal polygyny; easy tolerance of divorce and re-marriage and indeed 
minimal interference with sexual activities not constituting culturally 
defined “‘ incest”; some practice of the levirate and sororate ; approval 
of cross-generational marriages in both directions—i.e., older women 
also marry younger men with some frequency; tendency for repeated 
exchanges between biological or extended families in marriages (e. g., a 
man is more likely to marry the sister of his brother’s wife or cousin’s 


wife or the sister of his sister’s husband ). 


THE DATA 


Our calculations of inbreeding coefficients have been made from data 
summarized on a large genealogical chart showing the relationship and 
matings of 1118 members of the population during 8 generations 
(counting the initial, “ founding ” generation as one). The time period 
covered by the genealogies—as opposed to the present outbreeding popu- 
lation which has existed as such only since about 1870-1880—is from about 
1820 to September first 1948. This genealogical chart was compiled 
by Clyde Kluckhohn, Katherine Spencer, and J. N. Spuhler, and drawn 
by Ann Tulloch. (Unfortunately, the chart is too large to reproduce 
with this paper; copies are on deposit in the Peabody Museum, Harvard 
University, and the Institute of Human Biology, University of Michigan.) 
The population data basic to the Ramah Navaho census and the gen- 
ealogical chart were collected and worked up by members of the Ramah 
Project of Harvard University, especially Clyde Kluckhohn, Katherine 
Spencer, Janine Chappat Rosenzweig, and Nan Stoller. Individual pedi- 
gree charts (similar to that shown in figure 2) were traced by Helen 
Spuhler on a separate card for each inbred sibship in the population 
and these cards were used to calculate the coefficients for each inbred 
mating. 
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RAMAH NAVAHO MATING TYPES 


303 


Table 2 is a summary of the incidence of 62 different mating types 


observed in the Ramah population. 


Sixty-one of these mating types 


involve inbreeding, and of these 17 types involve matings where the 


parents are related in a single way, 44 types involve matings where the 


parents are related in two or more ways. 


Incidence of 62 mating types by sibship and siblings for 4 generations 


RELATION 
OF MATES 


], FIC 

9 HNeGA 
3. GUHNi 
4. FICIR 
5. HICIR 
6. F2C 

7. DF2C1IR 
§. H2C 

9. F2C1IR 
10. H2C1IR 
ll. F2C2R 
12. F3C 

13. H2C2R 
14. H3C 
15. F3C1IR 
16. H3C1R 
l7, H4C 


FIC 





Ramah Navaho Population. 


0625 
0312 
.0312 
0312 
0156 


0156 


0156 


0078 
0078 
0039 
.0039 
.0039 
.0020 


.0020 


0020 
0010 


0005 


0977 


0703 


0632 


.0352 


0352 


19 


*Includes one child of uncertain paternity. 
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— — 


See text for details. 


Totals 

Ss I 

13) (14) 
l 6 
1 5 
l 7 
3 22 
l +) 
9 49 
3 17 
] ] 
9 34 
4 19 
l 4 
B° 5 
l 5 
5 18 
7 15 
2 13 
] 3 
l 8 
] ] 
l 2 
l l 
2 3 


% of 

pop. to 
(15) 
5.076 
2.112 
2.956 
9.291 
1.900 

10.347 
3.590 
0.106 
3.590 
1.003 
0.211 
1.320 
0.135 
0.487 
0.406 
0.176 
0.020 


10.580 


0.952 


0.476 


1.429 
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oF MATES 


(1) 


35. 


36. 


42. 


ro 


46. 


(2) 
F2C 
F3C 
F3C1IR 


F2C 
DF3C 
H3C 
F3CIR 
H4C 
H4CIR 


F2C 


7. H2C 


H3C 


F2C 
H2C 


38. H3C 


F3C1R 
FC4 


. H2C 


F3C 


. H2C 


DF3C 


H2C 
F3C 
F3C1R 
H3C1R 
F4C 
F4C1IR 
H4C 
H2C 
F2C1R 
H3C 
F4C 
H2C 


. F3C1R 


H4C 
H4C1R 


. F2C1R 


F3C 


F2C1R 


5. F3C 


H3C 


F2C1R 
H3C 
H3C1R 


.0281 


0254 


.0283 


0117 


.0156 


.0166 


0186 


0105 


.0117 


.0137 


.0107 
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(9) (10) (11) (12) (18) 


3 


Totals 


I 
(14) 


9 


to 


wo 


% of 


pop. fy 


(15) 


2.619 


1.141 


1.375 


1.109 


0.211 


0.674 


0.252 


0.142 


0.145 
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TABLE 2—Continued 





RELATION GENERATION % a 
OF MATES p io 5 6 7 8 Totals Pop. i, 
S I S I S I S I S I 

(1) (2) (8) (4) (5) (6) (7) (8) (9) (10) (11) (12) (18) (14) (15) 

47. F2CIR 7 .0098 1 3 0308 
F3C1R 9 

48. F2C1IR 7 .0083 2 7 0.786 
H4C 1] 
F2C2R 8 

49. H2CIR 8 .0100 ] 3 0.406 
F3C1R 9 
H4CIR 12 

50. H2C1R 8 .0059 7 19 1.517 
F3C1R 9 

51. H2C1IR 8 .0049 ] 2 0.133 
H3C1R 10 

52. F3C 8 .0044 ] l 0.060 
H4C 1] 

53. F3C 8 .0042 2 ] 2 0.114 
H4C1IR 12 

54. H3C 9 .0029 l l 0.039 
H3CI1R 10 

55. H3C 9 0029 2 5 0.196 
F4C 10 

56. F3CIR 9 .0022 2 l 2 0.060 
H4C1IR 12 

57. H3C1R 10 .0015 l 6 0.122 
H4C 1] 
H2C 7 

58. DF3CIR 8 .0127 ] l 0.172 
H3CIR 10 
H2C 7 

59. H3C 9 .0107 l ] 0.145 
F4C 10 
H2C 7 

60. F3C 8 .0146 2 2 3 0.593 
H3C 9 
F4C 10 
HICIR 6 

61. F2C 6 .0342 2 l 2 0.926 
F3CIR 9 
F4C 10 

Totals, inbred matings 123 423 100.000 

62. “ Unrelated ” 0 155 496 vo" 











Totals, all matings 60 272 114 414 97 219 7 #14 278 919 





Pop. i, 


(15) 
0.398 


0.786 


0.406 


0.172 
0.145 


0.593 


100.000 
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In column 1 of table 2 a serial number is assigned to the various 
mating types. Column 2 gives the “ relation of mates” for each type. 


The following abbreviations are used in this column: 


A Aunt Ni — Niece 

C — Cousin R — Removed 

D Double U — Uncle 

F — Full 1 — First (or, preceding R, Once) 

G — Great 2 — Second (or, preceding R, Twice) 
H — Half 3 — Third 
Ne — Nephew 4 — Fourth 


For example, “ DH4C1R” means “ Double Half Fourth Cousins Once 
Removed.” 

Column 3 gives the power (p) of the exponent fixing the contribu- 
tion of each type of parental relationship to the inbreeding coefficient. 
For example, F1C (Full First Cousins) contribute (1/2)* to the coeffi- 
cient (as explained above). 

Column 4 gives the contribution for autosomal loci of each type of 
mating to the coefficient. This is the sum of (p) over each of the 
various sorts of relationship (or, over each inbreeding loop) between 
the mates involved in each mating. For example, in mating type 18, 
i= 0.0977, which is the sum of (1/2)*-+ (1/2)°+ (1/2)*. 

Columns (5-12) in table 2 give the incidence of inbreeding by sib- 
ships (S) and by individuals (I) in generations 5, 6, 7, and 8 of the 
Ramah population together with the totals (columns 13-14) for each 
mating type, and the proportion (column 15) which each inbred mating 
type contributes to the population mean inbreeding coefficient. 

Table 2, then, gives the distribution of all inbred individuals by 
mating type and generation. There are 123 inbred sibships containing 
423 individuals. No inbreeding is known for generations 1 through 4 
of the population, although the genealogies of the second, third, and 
fourth generations are established in fine detail. There are a total of 
38(+-) sibships with 199 offspring in generations 1-4. The number of 
sibships involved in the case of the 14 members of generation 1 is not 
exactly established. 

Included in generations 5-8 are 496 “unrelated ” individuals from 
155 “outbred ” sibships. There are, thus, a total of 278 matings of 62 
different types in generations 5-8. These matings produced a total of 


919 offspring. 
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Table 3 is a summary by generation of the inbreeding coefficients 
of the Ramah population, with information on all generations since the 
founding of the population. 

Some particular value of the inbreeding coefficient may have at least 
4 different referents (Ludwig, 44; Wright, 51). It may refer to: 
(a) an individual, (b) the offspring of a given mating, i.e., a sibship, 
(c) a generation, or (d) a breeding population. Table 2 gives the raw 
data for determination of the first three of these sorts of referents. 
Table 3 summarizes this data in different form and gives information 
on outbred matings of generations 1-8. 

Information is presented on 316 matings which produced 1118 off- 
spring in 8 generations. 193 of these matings were outbred and pro- 
duced 681 offspring. The remaining 123 matings were inbred and 
produced 423 offspring. 

The range of the inbreeding coefficient for individuals (or full sibs) 
is .0010 to .0977. The mean coefficient for all inbred individuals for 
generations 5 through 8 is .0175. The mean coefficient for the popu- 
lation (generations 5 through 8 only) is .0080. The mean coefficient 
for the entire population since the foundation (generations 1 through 8) 


is .0O66. 
GRAPHICAL REPRESENTATION OF THE POPULATION 


Homogeneous coordinates may be used to represent a Mendelian popu- 
lation (de Finetti, ’27, Haldane and Moshinsky, ’39). The frequencies 
of n genotypes which occupy a single autosomal locus may be represented 
by a point within a regular simplex in (n—1) dimensional space 
where the perpendiculars of its n bounding planes are coordinates whose 
sum is unity. Thus the population frequencies of the three genotypes 


(x AA + y Aa + zaa—1) formed by a single pair of autosomal alleles, 
A and a, may be represented (fig. 3) by a point P in an equilateral 
triangle Y)'Z, where z, y, z are the perpendiculars from P on YZ, XZ, 


and YY respectively. If we let the gene frequencies of A = p, and of 


a= (1—p) —q, then p—2—1/2y, and the points representing a 
population with a given gene ratio (p:q) lie on the straight line PQ, 





where @ is the projection of P on XZ. The equations for QZ = p are 


1/2y—p, and 2qz + (q— p)y — 2pz = 0. 
Under a system of random mating, z = p*, y = 2pq, and z = q’, and 
the points generated by a parabola y* — 4rz —0 (with vertex at the 


mid-point of the altitude of the triangle and with tangents to YY and 
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YZ at X and Z respectively) represent the genotype frequencies for 
autosomal single alleles in all possible ideally large random-mating 
populations. 

Under a regular system of inbreeding the frequency of heterozygotes 
is reduced by a factor (1—i) where i is, of course, the inbreeding 
coefficient. Thus the frequencies of the three genotypes in regularly- 
inbred populations are r = p* + ipq, y = 2(1 — i)pq, and z = q* + ipg. 


di 





2 





Q x 


Fic. 3. GRAPHICAL REPRESENTATION OF GENE FREQUENCIES AND GENOTYPE 
FREQUENCIES FOR THREE VALUES OF THE MEAN INBREEDING COEFFI- 
CIENT OF A POPULATION. THE PARABOLAS ARE DEFINED IN THE TEXT 

AND THE THREE Ponts P, P’, P” ARE DEFINED IN TABLE 4 
Three parabolas are shown in figure 3. The uppermost parabola 
represents the locus of points whose perpendicular coordinates measure 
the genotype frequencies for autosomal single alleles in a random-mated 
population. The middle parabola represents the locus of such points 
for the ideal case of a population practicing a regular system of in- 
breeding of such a degree that t= .0977. As shown above, this value 
is the maximum inbreeding coefficient observed for an individual Ramah 
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Navaho. The lowermost parabola represents the locus of points for a 
population following a regular system of full brother-sister mating for 
one generation. 

The parabola representing the locus of points for a system of in- 
breeding where i = .0080 (the mean for generations 5 through 8 of the 
Zamah group) would lie below but very close to the parabola for a 
system of random mating. It follows that the distribution of auto- 
somal genotype frequencies in the Ramah population would not be 
expected to show marked deviations from the proportion of heterozygotes 
expected under a system of random mating. We can illustrate this fact 
by considering the genotype frequencies for some particular gene fre- 
quency. Since the B gene (a*) of the ABO blood group system is absent 
in the Ramah population, we may treat the A(A) and O(a) genes as a 
pair of single autosomal alleles. (The distribution of the MN blood 
types at Ramah will be considered in the following section.) The popu- 
lation frequencies of the AO genes based on a pooled sample of 458 
individuals (Boyd and Boyd, 49; Mourant, ’50) are p= A = .126 and 
r= i = 874 (q =a? —0). 


TABLE 4 


Genotype frequencies for 4 ideal systems of mating and the observed 
gene frequencies of ABO blood group genes of A anda 
in the Ramah Navaho population 





Mating Inbreeding Group A Group O 

system Coefficient, i, x AA yAa zaa Total 
P 0 O16 220 .764 1.000 
Pos 0.0066 O17 .218 -765 1.000 
af 0.0977 027 .199 774 1.000 
td 0.25 043 165 .792 1.000 


The segment of the base of the triangle QZ represents the gene fre- 
quency p and the segment QX represents the gene frequency (1 — p) =r. 
The genotype frequencies in the population are determined jointly by 
these gene frequencies and the system of mating. We will consider four 
such systems. Three different mating systems are shown in figure 3 by 
points P, P’ and P”. The fourth mating system is shown only in table 4 
by point Pry and indicates the appropriate genotype frequencies for a 
system of mating where the mean inbreeding coefficient is i = .0066. 
Inspection of the values in the graph and in the table indicates that the 
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observed degree of inbreeding in the Ramah population would not be 
expected to bring about large deviations from the genotypical propor- 
tions theoretically appropriate for a system of random mating in the 
population. 


INDIRECT ESTIMATE OF THE INBREEDING COEFFICIENT 
OF THE RAMAH POPULATION 


Now that we have direct estimates of the inbreeding coefficient of 
the Ramah Navaho population, it is possible to compare the results with 
an estimate obtained by an indirect method. Haldane (’38) has shown 
that departures in the amount of heterozygosity from that expected 
under random mating can be interpreted to give evidence for the mating 
system in natural populations. 

Consider a pair of autosomal alleles, such as the genes responsible 
for the MN blood types in man, where heterozygotes are distinguishable 
by phenotypical inspection. We will represent the three phenotypes by 
M, MN, and N, and their respective genotypes by MM, Mm, and mm. 
Let the observed phenotypical frequencies in the population be 


aM+bMN+cN=n. 


If the population was reproduced by random mating (with regard to 
the alleles M and m) and the gametes from the previous generation are 
represented by the proportion u M:1m, the expected distribution of 


genotypes would be 
nu*/(u + 1)? MM: 2nu/(u +1)? Mm: n/(u + 1)? mm. 
The logarithm of the likelihood of the observation is 


L = (2a + b)log u— 2n log(u + 1). 
Thus 
dL /du = (2a + b)/u—2n/(u+1) =0, 
and 


u= (2a+ b)/(b+ 2c) + V2n(2a+ b)/(bd - o)*. 


If this value of u is applied to the expected distribution of genotypes 
given above, the expected and observed numbers for each phenotype may 
be compared. A x" test is appropriate. 

If the differences between expected and observed values are too large 


to be consistent with a hypothesis of random mating, and the departure 
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is in the direction of an excess of homozygotes, inbreeding is one of 
several possible explanatory factors. (In the case of the MN blood 
groups, probably we can safely exclude the possible explanatory factor 
of assortative mating. ) 

Bernstein (°30) found that expected ratios of the genotypes in an 
inbred population (with inbreeding coefficient 1) are 


(u? + 1u)MM:2(1—1)u Mm: (w+ 1)mm. 


Since the genotypes can be distinguished by serological tests, if we let, 


as above, the observed frequencies of phenotypes be 


aM+bMN+cN—n, 
then 
u== (2a+ b)/(b + 2c), 
and 
t= (4ac — b*)/(2a+ b)(b + 2c). 


The observed distribution of the MN blood types in a sample of 458 
Ramah Navaho (Boyd and Boyd, 49; Mourant, 50) is 


380 M+ 74MN+4N = 458 
u = 10.1707, the expected distribution is 
379.67 M +- 74.66 MN +- 3.67 N — 458.00, 
2, = 0.0358, and 0.8 << P < 0.9. 


Although the departure of the observed distribution from the expected 
values is in the direction of an increase in the proportion of homozygotes, 
that is, in the direction expected for an inbred population, the observed 
values are fully consistent with the hypothesis of random mating. But, 
as shown in the preceding section, the observed mean value of the 
inbreeding coefficient for the Ramah population is sufficiently small, 
and has been in effect for a sufficiently short period of time, that large 
deviations from the theoretical proportions appropriate for random 
mating are not to be expected. With this in mind, and noting that 
gene m is relatively rare in the population and that estimates of its 
frequency are subject to considerable fluctuation in small samples, it 
remains of interest to solve for the inbreeding coefficient. We find 
i= .0088. This indirectly obtained value is remarkably close to the 
directly obtained value .0080 for the mean inbreeding coefficient for al] 


sibs in generations 5 through 8. 
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DISCUSSION 


In previously published estimates of inbreeding coefficients in human 
populations, samples of matings have been classified into inbred versus 
outbred and a coetlicient of inbreeding estimated from the ratio of first 
cousin, second cousin, and other consanguineous matings (usually not 
beyond third cousins once removed) to “ unrelated’ matings. ‘The size 
of the sibship produced by the two sorts of mating is not entered into 
such estimates. We have 123 cases of 61 different sorts of inbred 
matings out of a total of 316 matings (including the sixty-second sort, 
“ unrelated”) for generations 5 through 8. 1% for inbred matings is 
(1.9847/123) — .0161 and for total matings 1 is (1.9847/316) — .0063, 
This estimate is lower than the estimate we obtained for total individuals 
in generations 5-8 (i9 = .0080) because our inbred families of genera- 
tions 5-8 have more offspring than the outbred families.’ 

Caution is required for the proper interpretation of our population 
mean coefficient of inbreeding. First, the value should be regarded as 
a minimum estimate of the true value of the coeflicient. As a result of 
extensive field work over a period of 16 years, the Ramah genealogies 
are remarkably complete for generations 2-8. There are some unknowns 
in the founding generation. If any of the 14 founders are themselves 
inbred, the mean values of the population coefficient would be slightly 
larger. However, since no inbreeding was practiced in generations 2-4 
(where genealogical information is “ complete”) this possible source of 
error would not greatly modify the total population estimates (on the 
assumption, of course, that there are no undetected relationships among 
the founders). 


+The crude data suggest that inbreeding in the Ramah population has not 
reduced fertility of inbred individuals, although when the data for all generations 
are pooled the average size of outbred sibships is 3.53 persons, and for inbred 
sibships, 3.44 persons. However, when the crude data are compared by generation, 
inbred sibships average consistently larger than outbred, the respective means for 
generations 5-8 being 3.44 as against 3.20. The larger sibships in the earlier 
generations (mean of generations 2-4 = 4.45) make a relatively large contribu- 
tion to the population mean. Of course such crude data may be highly mis- 
leading. Many of the families of generation 7 and 8 are incomplete, In a later 
paper on the demography of the Ramah yroup we will examine the size of inbred 
and outbred sibships in terms of exposure to pregnancy. For the present we may 
conclude that inbreeding in the present group has not been accompanied by 4 
marked decrease in fertility. 
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If there are undetected biological relationships among the 14 founders, 
the mean values for subsequent generations could be increased con- 
siderably. However, the known relationships among the founders is of 
the sort which might be expected on the basis of general knowledge of 
Navaho behavior. With due regard to the unknowns in the case, our 
opinion is that the estimated values for the population, while somewhat 
low, are not greatly low, perhaps not more than 10%. 

Another difficulty in arriving at an estimate of the true value of 
the inbreeding coefficient is the problem of determining true population 
size. Although the Ramah group is relatively isolated (see above) and 
relatively self-contained, there are a few individuals regarding whose 
membership in the population an arbitrary judgment must be made. 
Actually we made an independent classification of the individuals repre- 
sented in the census data of the Ramah area. These data include 
information of residential history, travel, participation in other Navaho 
communities, and work in extra-Navaho contexts. The census included 
information on a number of individuals who have a marginal partici- 
pation (in both biological and socio-cultural terms) in the Ramah popu- 
lation. Our result differed only with regard to a single outbred family 
with 9 children. JNS classed the entire family in the population while 
CK put the father in, and the remainder of the family out of the 
population. For the present calculations, all 11 individuals have been 
included in the population. If our criteria for membership were made 
extremely broad, the population size, as of September first 1948, would 
be increased by some 20 individuals, and the size over all generations 
reported here would be larger by about 70 individuals. Use of the larger 
population sizes would result in a slight reduction of the population 
estimates of the mean value of the inbreeding coefficients, but not a 
marked reduction because some of these 70 individuals are themselves 
inbred. 

In previously published estimates, data have not been available on 
the more remote relationships. It is of interest to compare our results 
based on “full” specification with those obtained when the coefficient 
is estimated from various degrees of incomplete specification. We will 
examine three cases: 1) estimates from the incidence of first cousin 
(1C) matings; 2) estimates based on the incidence of first and second 
cousin (1C + 2C) matings; and 3) estimates based (as in the Hiroshima 
and Kure data reported in Neel, et al., 49) on the incidence of matings 
involving first cousins, first cousins once removed, second cousins, second 
cousins once removed, third cousins, and third cousins once removed. 
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Percentage of total 
Estimate i, = (.0080) contributed 


Case Based on Relationships of iy by these matings 
5. FIC .0008 10.0 
2. FIC+F2C 0020 25.0 
3. F1C+ FICIR+F2C1R+F3C+ F3C1R .0034 42.5 


In the case of the Ramah population, our estimate for the generations 
reported here would have neglected about 57% of the known inbreeding 
had we restricted our study 1) to degrees of relationship equal to or 
greater than third cousins, and 2) to specification of a single (closest) 
type of relationship for each inbred mating. 

This datum is not of direct use in estimating the amount of under- 
estimation of the inbreeding coefficient in other studies. The datum is 
closely influenced by factors specific for the Ramah community, for 
example, the effective size of the local breeding population, and the local 
Navaho cultural patterns of mate selection. However, the datum is 
suggestive for the interpretation of reported values for the inbreeding 
coefficient. We might guess, for instance, that the estimate by Neel et al. 
(’49) and by Schull (753) in the Japanese case—where marriages of first 
cousins are culturally permissive and actually performed fairly often— 
would be closer to the true value, than the estimate of Haldane and 
Moshinsky (*39) based on data of Orel (32) in the case of Catholic 
marriages in the archdiocese of Vienna where marriage of first cousins 
is by special dispensation and not a culturally preferred type of mating. 
In these studies, the specification of relationships was restricted to 
degrees equal to or greater than third cousins and multiple relationships 
were not specified. 

SUMMARY 


Inbreeding coefficients are calculated for 316 matings which produced 
1118 offspring during 7 generations of the Ramah population. Of the 
316 sibships, 123 are inbred. The range of the inbreeding coefficient 
for individuals is .0010-.0977, the mean for inbred siblings is .0175, and 
the mean for the later generations 4 through 8 is .0080, and the mean 
for the total population, generations 1 through 8, is .0066. These values 
are minimum estimates for the population. 

Coefficients of the order observed for the Ramah Navaho, while about 
two times those found for Japan and 4 times those reported for Europe, 
are small compared to the estimate for the Dunker isolates in the United 
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States and compared to such regular systems of inbreeding as full first 
cousin mating. 

Although the small Ramah Navaho population is not highly inbred, 
mating of the more remote relatives makes a comparatively large con- 
tribution to the mean value of the inbreeding coefficient. 
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STUDIES ON THE VARIABILITY OF HUMAN 
URINARY EXCRETION PATTERNS 


BY H. ELDON SUTTON AND STEVEN G. VANDENBERG 
Institute of Human Biology, University of Michigan * 


INTRODUCTION 


S our knowledge of biochemical genetics increases, the question of 
A the genetic control of human differences comes more to the fore- 
front. Some of the “ inborn errors of metabolism,” such as alcaptonuria 
and phenylketonuria, have been recognized for many years. More 
recently attention has been focused on the differences which occur 
between people whose various metabolic processes fall within the normal 
range of variation. While some of these differences are undoubtedly a 
function of the environment of the individuals, evidence is accumulating 
which indicates that many of the observed differences are of a more 
fundamental nature and that some of them may have their origin in 
genetic differences between the individuals. 

Many experiments, although not designed to demonstrate these ge- 
netic relationships, show the marked individuality which exists for many 
body constituents. McGavack and Drekter (’45) found that thyrotoxic 
patients show changes in blood cholesterol after treatment with thiour- 
acil, but each person maintained his relative position with respect to the 
other patients. Groen and coworkers (52) subjected 60 individuals to 
diets containing from one to 900 mg cholesterol/day and were able to 
change the blood levels only slightly. They found that each person 
reacted in his own characteristic way to the different diets. They sug- 
gested that the cause of these individual peculiarities may be of a 
hereditary nature. 


* This study was supported in part by a grant from McGregor Fund of Detroit, 


Michigan. 
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Serum alkaline phosphatase in cattle has been shown to vary widely 
among individuals although remaining relatively constant for a particular 
animal (Crookshank, et al., 52); Wagner and Poindexter (’52) reported a 
constitutional individuality in the level of cholesterol-esterifying activity 
in human serum. Kajiro and coworkers (’52) have reported on a con- 
stitutional anomaly in which there is a lack of the enzyme catalase in 
the blood and certain tissues. Wilkins and Carlson (752), studying the 
qualitative differences in neutral 17-ketosteroids excretion in normal 
subjects, found that although the total amounts excreted in the urine 
might vary with environmental changes, the qualitative pattern of each 
person was highly characteristic for the individual. 

Among the more interesting recent studies in biochemical genetics 
have been the demonstrations that genetically controlled variations in 
human hemoglobin exist. Itano (’53) has summarized some of the 
properties of the different hemoglobins along with their relationships 
to certain blood conditions. 

Of the various anomalies of amino acid excretion, cystinuria has 
perhaps been most adequately studied. Dent and Harris (751) have 
summarized the characteristics of the different types of cystinuria, 
emphasizing the genetic relationships. The increased excretion by some 
people of the recently discovered amino acid B-aminoisobutyric acid is 
reported to have a genetic basis (Crumpler, et al., °51). Fink, Hen- 
derson, and Fink (751) also found the excretion of this amino acid to 
vary widely from person to person. With the advent of paper chroma- 
tography, a number of studies of individuality in urinary excretion 
patterns have been made (e. g. Williams, 51; Berry, Cain, and Rogers, 
‘51; Sutton, *51). An interpretation of these results in terms of the 
yenetic makeup of the individual has recently been made (Berry, *53). 

The results reported here are from a study of urinary excretion 
products of 32 individuals (16 sib pairs). The data were analyzed in 
such a way as to answer the following questions: (1) Does an individual 
on repeated sampling tend to show a characteristic level of excretion of 
the particular substances measured? (2) If so, does the individual tend 
to resemble his sibling in this respect more than would be expected by 
chance? (3) Can individuals be characterized by excretion patterns as 
opposed to excretions of particular substances? The affirmative answers 
which have been obtained should help form a sound basis for further 
study in the field of human biochemical genetics. 
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EXPERIMENTAL 


The subjects of this study were 32 patients of the Coldwater State 
Home and Training School, Coldwater, Michigan, and consisted of 16 
pairs of siblings. Two of the sib pairs, on the basis of their serological 
and morphological types, were probably identical twins (kindreds 08 
and 13). All of the patients were classified as mentally deficient with 
no known organic defect. Both members of kindred 29 were blind or 
nearly blind, however. There were 13 males and 19 females whose ages 
varied from 7 to 33 and whose 1Q’s varied from 28 to 83. The sexes 
are indicated in table 1. 

Three first morning urine samples were collected from each of the 
patients on alternate days. All samples were collected within a three- 
week period. The samples were frozen soon after collection and were 
stored in that state except when necessary to remove aliquots for analysis. 

Creatinine determinations were done by the method of Bonsnes and 

Taussky (°45) and inorganic phosphorus by the method of Fiske and 
SubbaRow (°25). All other substances were measured by paper chroma- 
tographic techniques using conditions previously specified (Sutton, ’51) 
except that aliquots of urine containing 40 yg creatinine were used for 
the bromcresol green (BCG) chromatograms. One alteration in the pro- 
cedure was the use of reflectance measurements to quantify those amino 
acids for which ninhydrin had been used for color development. For 
this purpose a Beckman model B spectrophotometer equipped with an 
integrating sphere reflectance attachment with beam expanding lens was 
used. The filter paper containing the ninhydrin-developed amino acid 
yas placed in the instrument so that the light beam (550 my) was 
centered on the point of maximum color density. Quantitative values 
for the unknowns were obtained by comparing their percentage reflec- 
tance with a standard curve based on the reflectances of known quantities 
treated in the same manner. Alanine was used as the standard for glu- 
tamine, ethanolamine and f-aminoisobutyric acid, and the values for 
these substances are expressed as mg alanine. This procedure requires 
considerably less time than most other methods in the literature and is 
capable of the same accuracy. 

The substances designated “ BCG acid .28,” “BCG acid .30,” and 
“ BCG basic .32 ” refer to spots of unknown composition which appear 


on bromcresol green chromatograms resolved in butanol-acetic acid-water 
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(80: 20:20) solvent.* The numerals indicate the Rf values at which 
these spots appear. “ Hippuric acid” and “ Lactic acid” refer to acid 
spots appearing on bromcresol green chromatograms at the expected 
position for the indicated acids. It is possible that other acids contribute 
to these spots, however. The substances designated “ DSA orange .85 ” 
and “DSA purple .90” appear on chromatograms which have been 
sprayed with diazotized sulfanilic acid. The spots occuring on the 
bromcresol green and the diazotized sulfanilic acid (DSA) chromato- 
grams have been quantified by carefully marking the area of the spots 
with a pencil and then measuring the area with a polar planimeter. 
In studies of the DSA spots, it has been found that over the range 
likely to be found in urine the area of the spot is approximately 
linearly related to the amount of material present. Occasionally the 
presence of salicylate metabolites prevented an accurate quantitative 
estimation of some of the substances detected chromatographically. In 
such cases, “ no value (NV) ” was recorded. 

The optical density measurements were made by comparing the 
undiluted sediment-free samples with water in a spectrophotometer at 
620 mu. 

RESULTS 


The original measurements of the urinary constituents are presented 
in table 1. The mean and standard deviation for each item are shown 
in table 2. In most cases, the measurements have been expressed as 
creatinine ratios in table 2, since this ratio is a more significant quantity. 
For most of these items, no population norms have been established ; 
however, no grossly abnormal samples were encountered and these mea- 
surements are comparable to those reported in other studies ( Berry, "53 ; 
Beerstecher, et al., 50; Young, et al., °51). The lack of a comparable 
control group makes it impossible to detect any minor peculiarities which 
may be associated with this experimental group. Studies on this group 
will, however, give information on individual consistency and sibship 
similarities. 

To test the reliability of the measurements as a function of individual 
differences, an analysis of variance was carried out considering the sample 
as composed of 32 individuals with three measures each. The F values 


The somewhat varying Rf values for these substances has resulted in dif- 
ferent designations by different authors. Berry (’53) designates as “ BCG acid 
30” what is here called “ BCG acid .28.” She does not report on the spot here 
called “ BCG acid .30.” 
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are shown in table 2 with their levels of significance indicated. The high 
level of significance obtained for nearly every item certainly justifies 
the conclusion that even items falling within the middle of the normal 
range occur at levels characteristic of the particular individual. 

A further hypothesis—whether or not siblings tend to resemble each 
was tested by considering 
the sample to be made up of 16 groups of two individuals each. For this 





other more than could be expected by chance 


purpose the averages of the three measures were used as the value for 
that individual. These results are also shown in table 2. The larger 
number of substances failing to show significance may in part be due to 
the fact that higher F values are necessary to attain a given level of 
significance when the sample is reduced in size as it is in the intersibship 
comparison. Further work with a larger number of subjects will be 
necessary before conclusions can be drawn regarding the intrasibship 
similarity. 

The origin of the intrasibship similarities would seem most probably 
to lie in their similar genetic makeup. The diet, which is the environ- 
mental variant most likely to influence the composition of the blood and 
urine, was fairly constant for all subjects because of their institutional 
status. Other studies (Thompson and Kirby, 749; Sutton, *51; Groen, 
et al., °52) have shown diet to be relatively unimportant in determining 
the levels of many body metabolites. The relative importance of heredity 
versus diet remains a subject for further study, however. Experiments 
are in progress in this laboratory which should throw additional light 
on this matter. 

In order to have a quantitative measure of the degree of intrasibship 
concordance, the indices shown in table 3 have been calculated. The 
similarity index is the reciprocal of the mean square average difference 
in standard scores between siblings for all the items measured. The 
column headed “r,.” is the product moment correlation based on the 
12 items of table 2 which showed an intrasibship similiarity significant 
at least at the 5% level. The column headed “ r.,” is the correlation 
coefficient based on all 23 items of table 2. Inspection of the table 
reveals a fairly close agreement between the two correlations and the 
index of similarity. Part of the difference is undoubtedly due to the 
fact that the correlation coefficient places emphasis on congruence of 
the extreme values while the index of similarity places equal emphasis 
on congruence of intermediate and extreme values. A choice between 


these indices is difficult until more is known concerning the relative 
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importance to heredity of absolute levels versus related patterns of 
excretion. It is interesting to note that neither of the two sets of twins 
(sibships 08 and 13) ranks highest in similarity as calculated in table 3. 

One interesting form in which individual variations occur is in the 
different patterns of urinary excretion. By this is meant the relative 
amount of each substance present when compared with other substances, 


TABLE 3 


Comparisons of profile similarities within sibships 


Index of 


Sibship Sex Similarity is rss 
23 eies 11.6 437 350 
22 2 9 10.6 563 163 
26 Jo 9.8 407 339 
10 2 g 8.9 432 374 
14 rote t 8.4 699 538 
13 2 9 8.1 625 501 
12 99 7. 132 042 
29 2 9 7.4 — .003 206 
25 niet 6.9 301 201 
11 a9 6.8 099 058 
09 J 3 6.8 396 236 
08 9 9 6.3 229 169 
28 ag 5.1 175 448 
07 9 9 oe 634 323 
27 99 3.3 180 - 308 
24 J 9 3.1 — .094 — .084 


Figure 1 shows such individual patterns or profiles. The values plotted 
are the standard scores, obtained for each item by subtracting the group 
mean from the individual value and then dividing by the standard 
deviation. The results thus expressed are comparable from one sub- 
stance to another. Sibships 22 and 23, whose profiles are shown in 
figure 1, possess a pattern observable to some extent in at least 10 of 
the 32 individuals studied. These similarities were first observed on 
graphs in which the items were plotted in a random sequence. In order 
to make the similarity more obvious, the individual items of the 4 
subjects whose patterns were most nearly consistent with each other 
(sibships 22 and 23) were ranked in order of decreasing size for each 
subject. The average rank of each item was then obtained, and the 
sequence of items shown in figure 1 is in terms of these increasing 
average ranks. 
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Fig. 1. Prorirces or URINARY EXCRETION PATTERNS OF Four Parrs OF SIBLINGS. 
THE SUBSTANCES ARE NUMBERED AS IN TABLE 2. 
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The resemblance between the profiles of sibships 22 and 23 is striking. 
By contrast, the profiles of sibship 14 seem to show a completely unrelated 
pattern, although the concordance within this sibship is high. The 
profiles of sibship 11 resemble those of sibships 22 and 23 for many of 
the items. The concordance of the siblings for the marked deviations 
which occur indicates that this deviation is not a matter of chance. At 
least three different patterns of excretion of the type discussed here have 
been found which were common to two or more sibships in this study. 


DISCUSSION 


Of particular interest in the past have been such anomalies of urinary 
excretion as alcaptonuria and cystinuria, conditions which are charac- 
terized by excretion of one or a few substances at levels many times 
higher than are usually found in a normal population. There probably 
are many such genetically controlled conditions whose presence goes 
unnoticed because the principal point of interference with metabolic 
processes lies too far removed from the easily measurable materials such 
as blood and urine. Many differences must also exist whose effect at 
any point in the metabolic processes is small. 

By the use of such an approach as is presented here, it is hoped that 
some inroad can be made into the area of the integral differences between 
individuals. The detection of differences in patterns of excretion should 
provide a new means of classifying human beings, perhaps into discrete 
groups, which will make family studies more fruitful. The primary 
variables which are responsible for these profile differences are as yet not 
well understood. A statistical approach such as multiple factor analysis 
should furnish considerably more insight in this area. 


SUMMARY 


The urinary excretion products of 16 sib-pairs who are mentally 
deficient were analyzed chemically for creatinine, phosphorus, amino 
acids, and other substances which can be detected on paper chromato- 
grams. Analysis of variance showed that most of the substances are 
excreted at levels characteristic of the individual and that siblings tend 
to resemble each other more than do unrelated persons. Several patterns 
of excretion were detected in the group. 
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MATHEMATICAL ANALYSIS OF THE GROWTH OF 
MAN, WITH SPECIAL REFERENCE 
TO FORMOSANS * 


BY SUNG-KEN QUO 
Department of Hygiene and Public Health, College of Medicine, 


National Taiwan University, Taipeh, Formosa 


INTRODUCTION 


HE growth of animals and plants consists of an increase of living 
it substance resulting in cell division. A human ovum which has 
been fertilized and implanted in the uterus takes about 280 days to 
d 
well known that the rate of human growth is greatest during the first 
year of life, and that the body of the human species, like that of other 


evelop into a baby which, in turn, continues to grow after birth. It is 


animals, ceases to grow after a certain age. The weight of a human at 
birth doubles itself in about 5 months and trebles itself by the end of 
the first year. 

The growth of the human body is three-dimensional, since it involves 
increase in the size of bones and other constituent organs. There has 
been a host of investigations on the growth of man which deal with the 
development of stature, chest girth, sitting height, span, weight and so 
forth. However, since the human body grows not only in one dimension, 
it is obvious that the growth of the body as a whole is most appropriately 
represented by the deevlopment of its weight, which is the resultant of 
the growth of the constituent organs. Therefore, body weight is the best 
index of growth. 

With regard to the development of weight in man at different ages, 
there have been a great number of studies, of which the results are 
summarized by Baldwin (’21) and Davenport (26A). Nevertheless, very 


+ Aided by a grant from The China Foundation, New York, N. Y. 











334 SUNG-KEN QUO 


few authors have attempted mathematically to analyze the development 
of body weight from birth to maturity. Although the process of human 
growth is of great complexity, it is conceivable that this process follows 
certain natural laws which can be expressed by mathematical formulas 
where the amount of growth y is the function of time z. Such formulas 
enable one to establish the norm of human growth, to foretell with 
considerable certainty the quantitative results of any future observations 
of a similar nature, and to form some idea of the causes of deviations 
between the calculated and observed results. 

Robertson (’15), as the result of certain theoretical considerations, 
found that growth in the weight of infants followed the autocatalytic 
equation. Lambolez (25) obtained an exponential equation to represent 
the law of growth in infants. Failing to analyze the developmental curve 
of body weight as it is, Davenport (26B) transformed it into a curve of 
annual increments of total body weight and concluded that this latter 
curve has two autocatalytic cycles, one being best fitted by a theoretical 
skew curve of Pearson’s type I and the other by a normal frequency 
curve. 

The curve of the development of weight from birth to maturity shows 
the same shape, roughly similar to a sigmoid, not only for Americans 
(Baldwin, ’21 and Watson and Lowrey, 52) and Nordics (Davenport, 
°26A and ’26B) but also for Chinese (Shirokogoroff, ’25) and Japanese 
(Yoshinaga, *30 and Department of Physical Education, ’27), except 
that the slope of the curve varies with race and sex. In view of the fact 
that no one so far has carried out any study on the physical growth of 
the Formosans from birth to maturity, it is of great interest first to 
collect as many data as possible for calculating the average weight of 
Formosans of both sexes at different ages, and then to analyze the weight 
curve to find the mathematical formulas that will serve as the norm of 
growth of the Formosans. Although the slope of the Formosan growth 
curve is different from that of others, the mathematical formulas obtained 
can be applied to other races when the various constants of the formulas 


are separately computed. 


PROBLEMS FOR SOLUTION 


In the mathematical analysis of the growth of Formosans, it is pro- 
posed that a number of problems which are important and interesting 
to various specialties of medical science and public health (physiology, 
anthropology, obstetrics, maternal and child health, pediatrics, school 
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health, internal medicine, and physical education) be solved at the same 
time. These problems and the specialties closely related to them are 
summarized in table 1. Although some of these problems may be readily 
solved if only the mathematical formulas for the entire growth curve 
are obtained, still further analysis of the curve is needed in order to 
find the solution to the other problems. 


TABLE 1 


Problems to be Solved 





a °° "PROBLEM 
I What course does the growth of man take in his lifetime? 
II What is the average birth weight of the Formosans? 
Iil What is the normal growth of Formosan infants? 
IV What is the normal growth of Formosan children of preschool 
age? 
V What is the normal growth of Formosan school children? 
VI What is the normal growth of Formosans from adolescence up 
to adulthood? 
Vil What is the maximum limit of growth of the Formosans? 
Vill At what age does the adolescence of the Formosans start? 
IX At what age does the rate of adolescent growth of the Formosans 


reach its maximum? 


xX At what age do the Formosans stop growing? 





MATERIAL 


Inasmuch as the physical growth of the Formosans has never been 
studied from birth to maturity, the weight data for the two sexes and 
different ages had to be collected from various sources which could be 
regarded as reliable. For the birth weight, the data obtained in 1942 
by Wu (°49) and those in 1950 by Chiu and his coworkers (751) have 
been utilized to calculate the mean. For the weight at each month of 
infancy, the data acquired by Wu (’49) have been used. For later years 
the data of the health survey collected in 1938 by the Department of 
Health of the Government-General of Formosa were used.? In all, data 


*It has been assumed that mean body weights of different persons at these 
various ages adequately represent the weight differences of a single group of 
persons successively observed at these ages. While some discrepancy might appear 
between these two types of observation, it would affect only the evaluation of the 
constants in the following formulas, and not the mode of analysis itself. 
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concerning 87,277 Formosan males and 83,376 Formosan females have 
been used in the present study. Table 2 presents the mean weight and 
the number of persons observed for each month of the first year and for 
each subsequent year of life. 

Figure 1 shows the growth curve of Formosan males and females 
plotted by using the average weights at different ages from table 2. 
Growth curves of Americans (Watson and Lowrey, ’52), Nordics (Daven- 
port, °26A and B) and Japanese (Yoshinaga, *30 and Department of 
Physical Education, ’27) are presented in figures 2 to 4 by the use of 
the data derived from the respective sources cited. The dots in all 
these curves represent the average weights actually observed. 


METHOD OF ANALYSIS 


On looking at the growth curves in figures 1 through 4, it is obvious 
that there are more than two cycles in the course of human growth 
from birth to maturity, although its curve looks roughly like a sigmoid. 
During the first year of life, there exists a cycle of rapid growth. The 
curve then begins to flatten out until about 5 years of age when it rises 
more rapidly again. It reaches a maximum of slope in this cycle at 
12 to 14 years and then gradually approaches the horizontal which 
indicates the state of maturity. 

The mathematical analysis of human growth consists in finding out 
from the observed results the mathematical formulas that express 
certain natural laws. If the human growth curve were a straight line, 
then the required formula would be y=a-+ bz, which is the general 
equation of a straight line; the values of the constants a and b can be 
found by simple calculation. However, as the growth of man demon- 
strates a complicated curve having more than two cycles, it is quite clear 
that a single formula is not enough to express it. There should be as 
many formulas as there are cycles involved in the whole curve. It is 
extremely difficult to find the formula for each cycle from the shape of 
the curve, since it is, as a rule, impossible to tell by mere inspection 
whether the portion of the curve belongs to a parabola, a hyperbola, a 
logarithmic curve, an exponential or a trigonometric curve. The required 
formula must in such cases be obtained by trial. Generally one has some 
sort of idea as to what kind of formula to expect, and in order to test 
the correctness of one’s expectations one can attempt to convert the 
curve into a straight line by means of certain artifices, such as trans- 


formation of the coordinates, use of a trial value for some constant of 
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AGE IN 


YEARS 
Birth 
1/12 
2/12 
3/12 
4/12 
5/12 
6/12 
7/12 
8/12 
9/12 
10/12 
11/12 


TABLE 2 


Observed Average Weight of Formosans 


Number of 
Persons 


Observed 


280 
100 
274 
317 
23 
519 
519 
576 
433 
337 
228 
190 
234 
648 
530 
459 
210 
033 
995 
914 
987 
S89 
798 
753 
545 
497 
157 
250 
299 
292 
124 
168 
973 
914 
819 
944 
338 
607 
6 704 


NNN NHONN Www Ww WwW 


NO =— = = — £9 1 PO PO PO PO OO PO 


8i 277 


MALE 
Average 
Weight in 


Kilograms 


3.131 

4.200 

5.250 

5.896 

6.462 

6.942 
263 
554 
835 
100 
318 
517 
718 
.003 
.769 
.194 
15.619 
16.931 
18.323 
19.808 
21.739 
23.561 
25.084 
27.608 
30.008 
33.409 
37.478 
41.389 
44.138 
46.905 
49.208 
50.884 
51.638 
52.444 
52.699 
52.830 
53.111 
53.318 
53.246 


—~oOoa@mmonn =! 


_ 
to 
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FEMALE 
Number of Average 
Persons Weight in 
Observed Kilograms 
215 3.065 
71 4.041 
264 4.995 
286 5.549 
455 6.014 
437 6.460 
476 6.853 
467 7.139 
288 7.427 
232 7.704 
140 7.973 
136 8.166 
149 8.438 
3 591 10.909 
3 376 12.203 
3 295 13.744 
3 124 15.000 
2 898 16.316 
2 925 17.869 
2 925 19.481 
2 850 21.334 
2 786 23.314 
2 656 25.489 
2 580 28.125 
2 523 31.474 
2 525 35.321 
2 239 39.690 
2 167 42.154 
2 321 44.036 
2 196 45.765 
2113 46.286 
2189 46.586 
1 905 47.220 
1 832 46.849 
1 843 47.138 
1 847 47.183 
7 971 47.719 
6 936 47.434 
6 147 


83 376 
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the formula, ete. If the attempt is successful, then all the constants of 
the formula assumed can be calculated by the method of least squares, 

This is but a first approximation of the constants to the theoretical 
curve. It is desirable, however, to obtain the best-fitting curve by 
successive approximations which will correct the computed values of all 
the constants and thus minimize the differences between the observed 
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points and the corresponding points on the fitted curve. If the second 
approximation to the best-fitting curve does not sufficiently minimize 
the deviations, we may proceed to the third and the fourth approxi- 
mations, and so forth, by repeating the same process, till the best fit to 
the observed curve is finally attained. 

It was found, by the procedure stated above, that the human growth 
curve has four cycles, of which the first and second cycles are best fitted 
by the catalytic equation 


IL 


y= k — e 


where y represents weight, x denotes time and k, a and b are constants 
to be determined. For the Formosans, the first evcle of growth ranges 
from birth to 9 months, and the second cycle from 9 months to 5 years. 

The third cycle of human growth, which was found to range from 
5 to 13 years in Formosans, is best represented by the compound interest 
equation 


y = eatbe 
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in which y and z denote weight and time respectively, and a and b are 
constants to be calculated. 

The fourth growth cycle of the Formosans extends from 13 years 
till maturity, at which the body ceases to grow. The best-fitting curve 
to the fourth cycle of human growth is the logistic equation 


where y represents weight, z denotes time, and d, k, a and b are constants 
to be determined. Among these four constants, d is the lower asymptote 
which denotes the total growth attained up to the end of the previous 
growth cycle, d+ k is the upper asymptote of the present cycle, b the 
inherent rate constant, and a the constant of integration dependent upon 
the position of the origin. The method of determining these constants 
is fully described in a previous paper of the writer (Quo, 750). 

The constants of the above theoretical equations obtained by the first 
approximation should be subjected to further approximations in order 
to get the best fit to the observed growth curve. The observed curve of 
a given growth cycle may be expressed by the equation 


y’ =f(u,v,w;2) 


where u, v and w are the unknown constants to be determined. Then 
the equation obtained from the first approximation is 


Y, = f(u, 01, 0,32) 


in which u,, v,; and w, are the first estimated values of the constants. 
Therefore 
UU +e; Vy +e; w—w44+e; 


where ¢,, ¢, and c, are the correction terms that will be calculated from 
a subsequent approximation. By the expansion of Taylor’s theorem, 
we get 
ss 
. of 0 
f(u, v,w;2) = f(u,,v;,w,; 2) + (u ui) = + (v—n) of 
: ‘ Ou, 


+ (w — w,)- ie. + 5[ w— u,)* ol ‘+ (v — v,)? or 


+ (w —w,)’ 2), +[ (w — v,)(w — w,) at 


o*f 
Ou,dw 


of 


+ (u— u,)(w — w,) — “hg (wu — u,)(v — 2) — 
7 1 1 


+ . 
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Neglecting terms of higher order than the first, we have 


,_y of _ of 4 of 
y =—_—=y— == C - Ce ~ ” Cs > 
Y , * Ou; Ov, Ow, 
which is a linear expression for the correction terms. Therefore, we fit 
regression 
y = yf — Y, = C2, + Cot + C2 
in which 


af af af 
y=, 9 to = - “it zr — '. “gy 
Ou; Ov; Cw, 


The values of c,, c. and c, are to be determined by the method of least 
squares. 
Now for the catalytic equation, we have 
f—k—es, 
therefore 
of 
ok, 
of , 
= oe —— 91 h? == —(k, —Y,), 
0a, 
of 
I, = - 
" Ob, 


where a,, b,, k, and Y, are respectively the estimated values of a, b, k 


=< rene = x(k, oor Y;), 


and y from the first approximation. 

Let Co =k. —k,; c, = a.—4,; Cc. = b, — by, in which k., a, and 
b, are the values of the constants that will be calculated from the second 
approximation. Then we fit regression 


y=y —Y, = Cy + C2, + CoX2 
where y’ denotes the observed weight. By the method of least squares, 
we choose Co, ¢, and c, to minimize 
f — =(y —— Co — 0,2, —— CaF 2 az 
Hence 
— 23) (y — co — ¢,2, — ¢22.) — 0, 
whence 


Co = Y — C,7, — Coke, 


in which g, 7, and Z, are respectively the means of y, z, and 2. 
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Similarly 3c -——=— 232, (y a Co Pe C12; Co) 0, 
1 
ee" 
and ——— 2>22(y¥ — Co — C12, — C2) = 0. 
( > 


On substitution for co, we get 
S21(y — Co — C2, — C2Z2) 
= Sril(y — 9 + ek, + Co€, — C12, — C2X2) 
= Y2i(y — 9) — 22. (4, — 41) — 2 Er, (t2 — Z2) 
=== (), 
Hence C22, (r, — £1) + ¢2B2,(22 — F2) = Ya,(y — 9) 
which is identical with 
¢; (2, — 4)? + eo (2, — #1)(@72 — #2) = S(21 — Z1)(y— FY) (A). 
Similarly 
Y22(y — Co — C2, — C2L2) 
= S2.(y — 9) — 4 S22(2, — £,) — 2D 72(22 — Fz) 
== (), 
Therefore ¢,Sr.(z, — %,) + ¢2¥22(2. — F2) = Sr2(y — 9) 


which is same as 
=) - - ' =) = \2 - - 
C12 (2; — Z,)(z, — Z2) ea Cod (a2 — Z2) = X(z2 a E2)(y ans ¥) (B). 
Equations (A) and (B) may be written 
C1811 > C2812 == Siy C1812 + C2822 = Say 
which are the normal equations whose solution gives the values of the 
correction terms c, and c,. The value of ¢c) can be easily derived from 
these terms. 
We may repeat this procedure, as well as those described below, for 
further approximations, till the best fit to the observed curve is obtained. 
For the second approximation to the compound interest equation, 


we put 


f— garda, 
of . 
Therefore t= = euthe o Y,, 
Ca, 
0 
Io= f =z reuthe = rY, 
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in which a,, b,, and Y, are respectively the calculated values of a, b and 
y from the first approximation. We fit regression 


y= y —Y,—c,2, + Cf, 


where c; = a, —4,;, Co =b, b,, and y’ is the observed weight. By 


the method of least squares, we choose c, and c. to minimize 


: . 
f= X(y — at, — c222)?. 
of 
Hence — == — 237, (y — 17, — C27.) — 0 
0c, 
of oy 
and = 2S22(y — ,2%,; — C.%2) = 0, 
Ol2 
whence ¢,D2%, + Sar, — Dey 
C1 2,02 + C227, = Sxry 
which may be written 
C1811 TT C2812 = Siy C1812 > C2822 = Say. 


These are the normal equations whereby the values of the correction 
terms c, and c, can be easily computed. 
For the second approximation to the logistic equation, we have 


f=<=d 
1+ e z 
hence 
—_ of l Y,—d 
os — ok, 1 4. p+ 1 k, , 
of key (Y,—d,)[k, —(Q 1,)] 

= - — 

5 Ca (1+ ¢ ) A 


le == ]. 


where Y,, d;, k,, a, and 5, represent respectively the estimated values 
of y, d, k, a and 6 from the first approximation. 
Let c, =k ky 3 Co == A a; ¢ b b, 3 ¢, = d, — d, in which 


ko, do, b. and d, are the values of the constants to be calculated from 


the second approximation. Then we fit regression 


, , i ' 
7 - Vy = €,2, + C2hq + Cyhg + C424, 


whe 
terl 
equ 
alre 


The 
mil 


cur 
pos 
a | 
Th 
sex 
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ind where y’ denotes the observed weight. The values of the four correction 


terms ¢,, C2, ¢; and c, can be calculated from the following normal 
equations which are obtained by the same method of least squares as 


already described : 

















By 
C1811 + C2812 + C3813 + C4814 = Sty 
C1812 + C2822 + C3823 + C4824 = Soy 
C1813 + Co823 + Ca833 + C4834 = Soy 
Ci814 + C2824 C3834 + C4844 == Say, 
These simultaneous equations are most readily solved by the deter- 
minants given below. 
Siy$12813814 $1181y513814 | $1181 281 y814 | 
SoyS228e3824 Si 282yS238e4 | S12S2282ySo4 
SsyS23833834 Si3S8sy833S34 | | SisS2sSsyS3z | 
] 
SsySo4834544 $148 4y834844 $1 4S824S8ayS44 
°.=- CC, = ———___ ?= . 
on $1181 2813814 | $11812513814 | $11812813814 
s SooSog8o4 81 28228e38o4 | | S; »Soo8 3824 
22823824 | 2382 282282382 
s s Se98a4 | 8; Sog8338 4 
$1 4824834544 | $1 4824834844 
| $318, 281381y 
Ss SeoSoe Soy 
S So S S 7, 
Si 48248345 ay 
C, = ——__—_—__. 
S Si2s S14 
i 
S S 2S So4 | 
Ss So S x | 
$1 4824834544 
RESULTS 
The mathematical analysis by the writer of different human weight 
es curves revealed that human growth, regardless of race and sex, is com- 
posed of 4 cycles, namely two catalytic cycles, a compound interest and 
h a logistic, passing from one cycle to the other without abrupt transition. 
m The general equation for each cycle is the same for all races and both 


sexes, but the constants of the equation vary due to the fact that the 


slope of the portion of the curve differs with sex and race. Slight 
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variations in the duration of growth cycles are observed among different 
races. 

The durations and the fina] equations for the 4 cycles of the 
theoretical growth curve of Formosans are shown below. These equa- 
tions enable one to predict the normal average weight for the two sexes 
at any age from birth to maturity. 


First Cycle: Catalytic Growth 
General Equation: 
Y= k — e2-be . 
J Fou 
or y = k — A e-be 
Formosans: 


Period: from birth to 9 months 


Equation : 





Male y = 8.8 — e}-7248-0.21722 (1) 
or y=8.8 — 5.6088¢-°- 21722 (2) 
Female y — 8.5 [oe e} 6870-0.1988z7 (3) 
or y= 8.5 — 5.4034¢-°-19882 (4) 
Second Cycle: Catalytic Growth 
General Equation: 
y — k et bz 
or y=k—Ae' 
Formosans : 
Period: from 9 months to 5 years 
Equation : Not 
Male y — 21.5 _—s e* 7390-0.1918¢ ( 5 ) 
or y = 21.5 — 15.4722¢-9: 19182 (6) 
Female y — 90.5 —— e* 6873-0.1958z% (7 ) 
or y= 20.5 — 14.6928¢-0-19882 (8) ' 
Third Cycle: Compound Interest Growth be o 
. . . aver 
General Equation: “4 
For 
y= garbe 
y pute 
or y= Ae (16 
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ent Formosans : 


Period: from 5 years to 13 years 


the : 
Equation : 
1a- 
Xeg Male y= e2-34034+0.0813852 (9) 
or y = 10.3851e° 8% (10) 
Female g am gh ee (11) 
or y = 9.4126¢°-4% (12) 
Fourth Cycle: Logistic Growth 
General Equation : 
k 
y—d+ oe 
1+ ev" 
vo. 
or y=—¢( — — 
: +r 1 + Ae™ 
1) 
Formosans : 
s) . 
") Period: from 13 years onwards 
3 . ; 
) Equation: 
4) 34.5 
? : me C ' a c +t se 
Male y — 19 “7. 1 + e 21071-0.447262 (13) 
34.5 
or y=—19 + —— ———_—_- 14 
? T TH 385590 an 
, — 25.5 = 
Female y — 7 , +. e} 06866-0 594742 (15) 
wae 25.5 16 
~ om 1+ se ate 
Note: y body weight in kilograms 
) xy == age in months for equations (1) to (4) 
) x age 1n years for equations (5) to (12) 
r=age in years —12 for equations (13) to (16). 
. The normal average weight of the Formosans at any time of life can 


be obtained by solving the above equations. Table 3 presents the normal 
average weight, both actually observed and theoretically calculated, of 
Formosans at each age. In the light of the small deviations of the com- 
puted values from the observed, it is clear that equations (1) through 
(16) excellently fit the 4 growth cycles of Formosans. 
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TABLE 3 





Average Weight of Formosans as Observed and as Theoretically Calculated 





AGE IN 
YEARS 


1/12 
2/12 
3/12 
4/12 
5/12 
6/12 
7/12 
8/12 
9/12 
10/12 
11/12 


MALE WEIGHT IN KILOGRAMS 
Deviations 
Observed Calculated (Obs.-Cale.) Observed Calculated (Obs.-Cale.) 


FEMALE WEIGHT IN KILOGRAMS 


Deviations 





4.200 
5.250 
5.896 
6.462 
6.942 
7.263 
7.554 
7.835 

8.100 

8.318 

8.517 

8.718 
11.003 
12.769 
14.194 
15.619 
16.931 
18.323 
19.808 
21.739 
23.561 
25.084 
27.608 
30.008 
33.409 
37.478 
41.389 
44.138 
46.905 
49.208 
50.884 
51.638 
52.444 
52.699 
52.830 
53.111 
53.318 
53.246 


3.191 
4.286 
5.167 
5.877 
6.447 
6.907 
276 
O74 
813 
096 
308 
518 
722 
10.947 
12.784 
14.302 
15.597 
16.919 
18.353 
19.908 
21.595 
23.426 
25.411 
27.565 
29.967 
33.545 
37.379 
41.103 
44.394 
47.065 
49.089 
50.543 
51.549 


aoeao onsen 3 


2° 


52.675 
52.968 
53.359 
53.485 
53.498 





—0.060 
0.086 
0.083 
0.019 
0.015 
0.035 
0.013 

—0).020 
0.022 
0.004 
0.010 
0.001 
0.004 
0.056 
0.015 

—0.108 
0.022 
0.012 

—0.030 
0.100 
0.144 
0.135 
0.327 
0.043 
0.041 
0.136 
0.099 
0.286 
0.256 
0.160 
0.119 
0.341 
0.089 
0.217 
0.024 

—0.138 
0.248 
0.167 


0.252 


3.065 
4.041 
4.995 
5.549 
6.014 
6.460 
6.853 
.139 
427 
704 
973 

8.166 

8.438 
10.909 
12.203 
13.744 
15.000 
16.316 


J J -3 =] 


3.097 
4.071 
4.870 
5.524 
6.061 
6.501 
3.861 
157 
399 
597 
O19 
223 
421 

10.569 
12.335 
13.788 
14.981 

16.293 
17.854 
19.564 
21.437 
3.490 
5.740 
8.205 
1.784 
5.519 
39.127 
42.083 
44.197 
45.565 
46.395 
46.878 
47.153 
47.307 
47.393 
47.44] 

47.490 
47.499 
47.500 


oo 


Soon s =J 


2 2 


0.032 
-0.030 
0.125 
0.025 
—0.047 
—0.041 
—).008 
0.018 
0.028 
0.107 
0.046 
0.057 
0.017 
0.340 
0.132 
0.044 
0.019 
0.023 
0.015 
0.083 
0.103 
0.176 
0.251 
0.080 
0.310 
0.198 
0.563 
0.071 
0.161 
0.200 
0.109 
0).292 
0.067 
0.458 
0.255 
—).258 
0.229 
—0.065 
—0.216 
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DISCUSSION 

Growth consists of an increase in the number of tissue cells as a 
result of cell division. Hence it follows that growth takes place according 
to geometric progression in some form or other. Food material, ingested 
and assimilated in the body, is absorbed into the blood and carried to 
the various tissues of the body where it is partly used up, and partly 
transformed into the components of the tissue cells, such as nuclear 
material and cytoplasm, which are utilized in cell division to build up 
new cells. This biochemical transformation of food material is caused 
by the catalytic action of enzymes, which have vitamines as active com- 
ponents and are regulated by hormones. 

; 

chemical reaction by catalysis, the reaction velocity = at any 


instant is proportional to the amount of material present at that 


ly 


moment, i. e., 
dy 


dz vin b(k—y) 


where y represents the amount of material transformed at time z, k 
denotes the original amount of material, b is the proportionality constant 
and / — y equals the amount of material which remains untransformed 


at time x. This differential equation may be written 


dy 
: == bdz. 
k—vy 
Integrating, 
dy 7 
= b dx 
J i y . 
we get log, ( I: ~ y) + C= bz + c’, 
or log, (i - -y)= br + c—C’, 


where c and c’ are constants of integration. 


Put a een 
then log-(k — y) =a— br. 

Hence y= I; — et-be. 

= y = k — Ae (A = e*), 


both equations expressing the two catalytic cycles of human growth. 
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Since growth is due to geometric increase in the number of tissue 
cells, it can also be considered from another aspect of the chemical 
reaction which occurs in the tissue to synthesize nuclear material and 
cytoplasm. The reaction velocity sy at any instant z is also propor- 
tional to the amount of the product, y, formed by the reaction at that 
moment. If this proportionality is converted into an equality by inser- 
tion of a proportionality constant b, we obtain the differential equation 
for the rate of formation of the product by the reaction, namely, 


dy 
“ == by. 
dz 
Separating the variables, 
dy 
~ == bdr, 
y 
and integrating, 
* dy F 
f — = h dr. 
e y . 
we get log-ytec=br+c’. 
Putting a=c —c, 
we have log. y=a-t be. 
Therefore y =m etrbz, 
or y= Ae (A == ¢%). 


which are the equations expressing the compound interest law in Nature. 
The third cycle of human growth follows this law. 

The logistic equation, discovered by Verhulst (1838, 1845, 1847) and 
independently by Pearl and Reed (°20, °22, °23), gives an S-shaped 
curve with an upper and a lower asymptote and describes accurately 
the growth trend of human populations which cannot surpass its upper 
limit in a growth cycle. The fourth cycle of human growth is best 
expressed by the logistic equation. It is of considerable interest that 
the growth of human individuals frog adolescence to maturity follows 
the same logistic law as the growth of human populations. 

The maximum limit to the growth in man is represented by the 
upper asymptote of the logistic curve, d +k, of which d is the lower 
asymptote of the fourth cycle and & corresponds to the total amount 
of growth attained in the same cycle. Therefore, the maximum limit 
to the growth of the Formosans for each sex is 
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Male d+k 19+ 34.5 = 53.5 Kilograms, 
Female d- k= 22+ 25.5 = 47.5 Kilograms. 


Since the logistic growth cycle extends from adolescence to adult- 
hood, we can calculate the age of maximum rate of adolescent growth 
by finding the coordinates of the point where the rate of growth is a 
maximum. ‘This point is the point of inflection of the logistic curve. 
After that point is passed, the rate of growth declines progressively, till 
finally the curve stretches along horizontally in close approach to the 
upper asymptote, over which the body cannot grow. Now the logistic 
equation to the fourth cycle of human growth is 

k 


ym dt ete 


We put g—y—d=— . er 
lata 
dx (1+ e%#)? 
kb — kb — kheothe 
Bie (1 — error)? = 
bz? — kbz 
-— 


then 


which is the equation to the first derivative curve. Therefore 
sf 1 , 
2” = ; (2bz — kb) 2 —0. 


, 


At the point of inflection evidently 2’ is not equal to zero, hence 


k 
2 = ~e 
k 

and y=s+ d. 
Therefore 1 + etre a 2. 

1 e 
anc .-_ — " 

b 


Accordingly, the age at which the maximum rate of adolescent growth 


occurs in the Formosan male is 
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a 1.21071 
———— © om | 2 — 4 x 
z+12=— b + 12 0.44726 i 12 = 14.707 


= 14 years and 8.5 months. 


Similarly, the age of maximum rate of adolescent growth for the 
Formosan female is 

a 1.06866 

xz 12 = — 12 = —— 

+ b Tr 0.59474 


= 13 years and 10 months. 


+ 12 == 13.797 


With a view to comparing the growth velocities at different ages, 
the curves of annual increments in the body weight of the Formosans 
are illustrated in figure 5. 

After birth the velocity of growth begins to diminish at first rapidly 
and then gradually until between 4 and 5 years, both sexes always 
running roughly parallel with each other. It then rises by degrees, 
in such a fashion that the female exceeds the male all the time. At 
between 11 and 12 years in the female and between 12 and 13 years in 
the male, the velocity of growth rises abruptly and reaches a maximum 
at between 13 and 14 years in the female and between 14 and 15 years 
in the male, though the latter maximum is higher than the former. 
The point of maximum growth velocity corresponds to the point of 
inflection of the logistic cycle. It is conceivable that the starting point 
of this abrupt rise in growth velocity indicates the beginning of adoles- 
cence. The age of maximum rate of adolescent growth for each sex 
which has been theoretically calculated, can also be estimated graphically. 

After this peak is reached, the velocity of growth diminishes at first 
slowly, then rapidly, and again slowly, till it practically reaches zero 
which indicates maturity or adulthood. After the age of maximum 
rate of adolescent growth is passed, the males always weigh more than 
the females. It is obvious from figure 5 that the age at which growth 
stops is 26 years for Formosan females and 31 years for Formosan 
males. Further gain in weight after these years may be considered as 
due solely to the deposition of fat, which does not involve increase in 
the number of tissue cells resulting from cell division. The curves in 
figure 5 make it clear that females mature earlier than males. 

It is well known that the endocrine glands play a vital role in 
growth by a specific effect of their hormones on growing tissue and by 
their influence on the body’s metabolism. The anterior pituitary 
secretes a specific growth hormone whose principal actions are accelera- 
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tion of bone growth and tissue building. Less directly the gland 
influences growth by the action of its hormones on other glands of 
internal secretion. The thyrotropic hormone of the anterior pituitary 
regulates the secretion of thyroxine of the thyroid gland, which is prob- 
ably second only to the pituitary in relative importance in its influence 
on growth and development. A normally functioning thyroid pro- 
motes skeletal growth, increases the metabolic rate and promotes sexual 
maturation. 

The gonadotropic hormones secreted by the anterior pituitary cause 
the ovarian cycle and the secretion of ovarian hormones in the female, 
and bring about the spermatogenesis as well as the secretion of testis 
hormone in the male. The gonads are also necessary for proper develop- 
ment of the body. The importance of the parathyroid glands to normal 
growth and development can be appreciated from the fact that these 
glands are responsible for normal bone development. 

It is well established (DeLee and Greenhill, ’48) that the functions 
of the maternal anterior pituitary and thyroid glands are increased 
during pregnancy and gradually return to normal after childbirth. The 
fetus which has been growing at an accelerated rate under the influence 
of the various growth-inciting hormones from the mother’s blood, slows 
down the velocity of its growth after birth (see figure 5), presumably 
for the reason that these hormones come to the infant only insufficiently 
from the milk. As time goes by, the child’s own endocrine glands 
gradually increase their functions, so that the velocity of growth rises 
correspondingly. 

When the function of the anterior pituitary and thyroid glands is 
increased so greatly as to stimulate the sex glands, the velocity of growth 
is sharply accelerated. This is the onset of adolescence which is to be 
followed before long by puberty or sexual maturation. As can be seen 
from figure 5, the age of onset of adolescence is between 11 and 12 years 
in the Formosan female and between 12 and 13 years in the Formosan 
male. 

Puberty begins when the sex glands become mature as a result of 
continuous stimulation by the anterior pituitary and thyroid hormones. 
The first menstruation is generally recognized as an objective indicator 
of puberty in the female, while in the male there is no marked sign, 
except for the onset of nocturnal emissions which does not necessarily 
coincide with the onset of puberty. According to S. R. Ko (’33) and 
K. E. Ko (’33), the average age of the first menstruation in the For- 
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NO. 


VIII 


PROBLEM 


TABLE 4 


Solution of the Problems Proposed 





First Cycle: Catalytic Growth 


Period for Formosans; from birth to 9 months 


SOLUTION 





Second Cycle: Catalytic Growth 


Period for Formosans; 


Period for Formosans; 


Fourth Cycle: Logistic Growth 


Period for Formosans; 


Birth Weight: 


Maximum Limit 
of Growth: 


Age of Onset of 
Adolescence : 


Age of Maximum 
Rate of 
Adolescent 
Growth: 


Age at which 
Growth stops: 


Male 


Female 


see 


Male 


Female 


Male 


Female 


Male 


Female 


Male 


Female 


from 9 months to 5 years 


Third Cycle: Compound Interest Growth 
from 5 years to 13 years 


from 13 years onwards 


3.191 Kilograms 


3.097 Kilograms 


table 3 


53.5 Kilograms 


47.5 Kilograms 


(see table 3) 


between 12 and 13 years 


(see figure 5) 


between 11 and 12 years 


14 years and 8.5 months 


(see figure 5) 


13 years and 10 months 


31 years 


26 years 





(see figure 5) 
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mosan female is 14 years to 14 years and 8 months. This is nearly the 
age of maximum rate of adolescent growth in the Formosan female that 
has been mentioned previously. It would seem that at the onset of 
puberty the velocity of the adolescent growth is at its maximum. In 
other words, the advent of pubescence is preceded by an acceleration, 
and followed by a retardation of growth in weight in the two sexes. 

The problems proposed in table 1 are already solved by the results 
and discussion in the present study; the solutions to these problems 
are summarized in table 4. 


SUMMARY 


In order to solve problems which are of practical importance in 
various branches of medical science, a mathematical analysis of the 
human growth in weight has been carried out with special reference to 
the growth of the Formosans. Data about 87,277 Formosan males and 
83,376 Formosan females, collected from 1938 to 1950, have been used 
in the present investigation. 

It has been found that the growth of man has 4 cycles from birth 
to maturity. Both the first and the second growth cycles obey the law 
of catalytic reaction, except for the difference in rate of growth between 
the two cycles. The third growth cycle is ruled by the law of com- 
pound interest in Nature, whilst the fourth growth cycle follows the 
logistic law that describes the growth of human populations. These 
4 cycles pass from one to the other with gradual transition. 

The mathematical procedure for the analysis of the human growth 
curve has been thoroughly described, and the periods and equations of 
the 4 cycles for the two sexes of the Formosans have been calculated 
with accuracy. 

The body weights at different ages of the Formosan male and female 
computed from the equations obtained were found to yield excellent 
agreement with the observed data from birth to maturity. 

The ages for the onset of adolescence and the maximum rate of 
adolescent growth among Formosans of both sexes, as well as for physical 
maturity, have been determined; and the maximum limits of growth 
that Formosan males and females may reach in their lifetime, have also 
been calculated. 

The changes in the growth velocity throughout the whole period of 
body development have been discussed from the endocrinological point 
of view. 
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